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Introduction 


The  1 9 8 3 /8 4-?r e s e a r ch  year  saw  the  completion  of  the  ONR 
funded  research  ab  UMIST>on  turbulent  convective  heat  transfer  in 
180®  bends.  As  last  year,  the  main  report  is  kept  very  brief  but 
technical  papers  and  reports  are  included  as  appendices.  The  main 
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research  effort  during^  the  year,  b-M  been  directed  at  bends  of 
circular  cross  sectloriT^,  X®ctlons  3  and  4  deal  respectively  with 
experimental  and  computational  aspects  of  this  work.  These  are 
preceded  Sectl-oir-2;>  by  a  brief  account  of  further  computational 
work  in  the  square-sectioned  duct .y,  Finally ,  Section  5  attempts  a 
summing  up  of  what  has  been  achieved  and  the  questions  that 
remain . 


In  Appendix  1  a  complete  list  of  theses  and  external 
publications  arising  from  this  project  is  provided.  Here  we  have 
included  contributions  made  jointly  with  Professor  Humphrey's  team 
at  UC  Berkeley  with  whom  the  whole  research  programme  has  been 
conducted  in  cordial  and,  we  believe,  effective  collaboration. 


2.  Computational  Studies  on  Square-Sectioned  180°  Bend 

Following  the  completion  of  Johnson's  [1]  PhD  work  Professor 
Choi,  an  academic  visitor  to  UMIST,  continued  explorations  into 
the  cause  of  the  serious  disagreement  between  experiment  and 
computation  for  the  case  of  a  square-sectioned  180°  bend.  Several 
modifications  were  Introduced  to  the  numerical/physical  model 
adopted  by  Johnson. 

(i)  Initially,  better  estimates  of  the  mean  velocity 
profiles  and  the  corresponding  turbulence  energy  and 
dissipation  profiles  at  the  station  just  upstream  of  the 
bend  were  attempted.  In  particular,  the  representation 
of  the  weak  streamwise  vorticity  generated  by  the  aniso¬ 
tropic  turbulent  stress  field  during  the  flow  develop- 
men  in  the  straight  30  diameter-long  entry  tangent  was 
considerably  Improved  and  this  led  to  significantly 
better  agreement  with  experiment  at  the  3°  position 
(l.e.  just  after  the  start  of  the  bend).  There  was. 
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however,  negligible  effect  of  these  changes  in  entry 
condition  at  stations  further  downstream.  The  resultant 
flow  pattern  at  90°  seemed  nearly  insensitive  to  what 
conditions  were  prescribed  three  diameters  upstream  of 
the  bend  inlet. 

(ii)  The  PSL  technique  12  3,  facilitating  a  very  fine  near¬ 
wall  mesh  (without  core  penalty)  was  introduced  first  on 
the  flat  top  wall  and  subsequently  on  the  curved  inner 
and  outer  walls  of  the  bend.  The  desirability  of  doing 
this  is  seen  in  Figure  1  which  shows  that  the  maximum 
secondary  velocity  occurs  very  close  to  the  top  wall. 
With  the  usual  coarse-grid  "wall  function"  treatment 
the  secondary  flow  is  appreciably  'clipped*. 

Significant  refinement  of  the  primary  mesh  was  also  made 
both  in  the  streamwlse  and  cross-stream  directions.  The 
finest  grid  employed  in  the  computations  was  27  x  47  x 
130.  Considering  that  the  third-order  QUICK  scheme  was 
used,  it  was  our  view  that  numerical  errors  had  been 
reduced  to  unimportant  levels. 

Computations  showed  that  (principally  due  to  the  use  of 
PSL)  the  predicted  behaviour  at  «5°  into  the  bend  was 
now  in  satisfactory  agreement  with  experiment.  Figure  2, 
if  allowance  were  made  for  what  appeared  to  be  an  error 
in  the  measured  bulk  flow  rate.  In  Johnson's  work, 
however,  the  main  disagreement  with  the  measured 
behaviour  had  not  arisen  until  90°  into  the  bend.  From 
Figure  3  it  is  seen  that  while  definite  improvements 
have  resulted  from  the  numerical  refinements  summarized 
above,  there  still  remains  very  serious  disagreement 
between  the  experimental  and  predicted  behaviour. 

(lii)  The  conclusion  drawn  was  that  shortcomings  in  the 
turbulence  model  were  principally  responsible  for 
remaining  discrepancies  between  the  computed  and 
measured  flow  pattern.  As  a  step  in  testing  this 
conclusion,  a  modification  was  introduced  to  the 

"standard"  form  of  the  dissipation  rate  equation.  With 
the  k ~ £  eddy  viscosity  model  (employed  over  the  whole 
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duct  except  the  PSL  region)  the  source  of  e  is  taken 
equal  to 
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where  \>t  is  the  kinematic  turbulent  viscosity  and 
denotes  the  mean  velocity  in  direction  x^.  In  some  of 

the  Pi's  earlier  studies  it  has  been  found  that  dropping 
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the  part  of  this  term  involving  J  brought  at 

least  qualitatively  the  correct  response  of  vt  to 
streamline  curvature.  Accordingly,  a  set  of  180°-bend 
computations  was  made  for  the  case  where  the  dissipation 
source  term  was  abbreviated  to 
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We  note  that  In  a  straight  thin  shear  this  form  still 
gives  essentially  the  same  behaviour  as  the  "standard" 
version  since  the  omitted  term  makes  a  negligible 
contribution.  For  the  strongly  curved  180°  bend, 
however,  this  model  variant  predicts  major  differences 
in  the  flow  pattern  (Figures  *l  and  5)  including,  at  90°, 
large  troughs  in  mean  velocity  near  the  inside  of  the 
bend  qualltati vely  in  line  with  experiment.  Overall  the 
computed  flow  pattern  cannot  be  said  to  be  in  any  better 
agreement  with  experiment  than  when  the  "standard"  model 
was  employed;  nevertheless,  the  results  did  at  least 
underline  the  sensitivity  of  the  predicted  flow  to  the 
turbulence  model  employed.  Birch  [3],  in  a  parallel 
study  employing  a  simple  algebraically  prescribed 
turbulence  model,  has  drawn  analogous  conclusions. 


In  continuing  work.  Professor  Choi  (now  returned  to  Korea 
University)  is  trying  to  obtain  predictions  to  this  flow  using  an 
algebraic  stress  model  (ASM),  a  level  of  turbulence  modelling 
which  (while  considerably  more  complex)  is  known  to  give 
more  reliable  predictions  in  complex  flows  than  models  based  on 
the  (Isotropic)  eddy  viscosity  concept. 
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3 .  Experimental  Work  -  Circular-Sectioned  180°  Bend 

During  the  present  contract  year  the  construction  and  shake- 
down  of  the  180°  bend  apparatus  has  been  completed,  the  heat-flux 
meters  (which  were  at  a  prototype  stage  at  the  time  of  the  last 
report)  brought  to  their  final  stage  of  design  development  and  a 
substantial  series  of  experiments  undertaken.  Manuscripts  have 
been  prepared  on  both  the  design  and  functioning  of  the  heat-flux 
meters  and  on  the  180°  bend  test  programme;  these  are  included  as 
Appendices  2  and  3  respectively. 

As  noted  in  last  year's  report,  it  did  not  appear  feasible  to 
provide  a  high-quality  experiment  with  a  uniform  heat  flux 
boundary  condition.  Accordingly,  a  uniform  temperature  rig  was 
designed.  Thus,  the  major  instrumentation  requirement  was  for  an 
accurate  heat-flux  meter.  The  instrument  which  provided  the  basis 
for  determining  local  values  of  Nusselt  number  in  the  U-bend  is 
shown  in  Figure  6.  A  set  of  five  such  gauges  set  at  45°  intervals 
around  a  section  of  straight  pipe  showed  a  reproducibility  for  any 
gauge  within  2%  in  successive  tests  and  a  maximum  variation  among 
the  five  gauges  of  5%  (which  is  probably  no  more  than  the  possible 
circumferential  variation  in  Mu  due  to  flow-field  non-uniform¬ 
ities)  . 

In  the  U-bend  itself  25  such  heat-flux  meters  were  installed 
(five  different  circumferential  angles  at  five  stations)  with  a 
further  five  at  a  (movable)  station  downstream  of  the  bend.  They 
have  proved  mechanically  robust  and  highly  accurate;  at  a  Reynolds 
number  of  1.1  x  10^  deduced  values  of  Nusselt  number  are  estimated 
to  be  accurate  within  ±5>.  A  mapping  of  both  the  heat-flux  field 
and  the  associated  temperature  field  within  the  air  in  the  tube  at 
the  same  stations  was  provided  at  a  Reynolds  number  of  about 
60,000.  This  is  close  to  that  for  which  Azzola  and  Humphrey  [4] 
have  obtained  laser-Doppler  velocity  measurements  in  an 
identically  proportioned  apparatus.  Figure  7  shows  the  circum¬ 
ferential  variation  of  Nusselt  number  measured  at  the  six  stations 
considered.  By  the  90°  station  (i.e.  halfway  around  the  U-bend)  a 
3:1  variation  in  Mu  has  developed  between  the  outside  and  inside 
of  the  bend.  There  are  two  sources  of  this  variation:  it  is  due 
partly  to  the  secondary  flow  depositing  hot  fluid  on  the  inside  of 


the  bend  and  partly  to  the  augmentation  of  turbulence  mixing  near 
the  outside  of  the  bend  and  its  damping  near  the  inside  that 
results  from  the  different  signs  of  the  angular  momentum  gradient 
in  the  two  cases.  What  Initially  seems  curious  is  that  six 
diameters  downstream  of  the  bend  there  is  still  a  2:1  variation  in 
Nu  despite  the  fact  that  the  streamlines  are  now  virtually 
straight  and  the  secondary  flow  is  very  small.  The  reason  turns 
out  to  be  that  the  hot  fluid  that  had  accumulated  on  the  inside  of 
the  bend  further  upstream  simply  stays  in  the  same  (circum¬ 
ferential)  position  depressing  the  rate  of  heat  transfer. 

In  further  work  undertaken  since  the  termination  of  the  grant 

the  corresponding  distributions  of  Nusselt  numbers  at  Reynolds 
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numbers  of  1.1  x  1 0  and  2.0  x  10  have  been  obtained;  these  data 
are  also  reported  in  Appendix  3.  Future  plans  call  for  additional 
experiments  with  thin  inlet  boundary  layers  and  cases  of  low 

4 

Reynolds  numbers  (<  10  )  where  laminarlzatlon  of  the  flow  on  the 
inside  of  the  bend  may  occur. 


M .  Computational  Work 

Last  year's  report  described  the  toroidal-coordinate  semi- 
elliptic  procedure  developed  for  analyzing  flows  in  circular- 
sectioned  bends.  At  that  time  it  had  been  applied  to  laminar 
flows  in  180°  bends,  ref.  C5],  and  turbulent  flow  in  90°  bends  [63 
employing,  in  the  latter  case,  the  standard  k ~ t  eddy-viscosity 
model  in  the  main  part  of  the  flow  and  a  mixing-length  model  in  a 
thin  annular  ring  adjacent  to  the  pipe  wall. 

Subsequently,  turbulent  flow  computations  using  precisely  the 
same  model  were  also  made  of  the  180°  bend  flow  being  studied 
experimentally  at  UC  Berkeley.  The  numerical  calculations  in  fact 
preceded  the  experimental  data,  so  when  the  latter  became 
available  there  was  some  satisfaction  at  the  encouraging  degree  of 
agreement.  The  manuscript  on  this  work  is  included  as  Appendix  <t. 
Figure  8,  which  is  drawn  from  this  paper,  shows  measured  and 
computed  profiles  of  the  mean  streamwise  (primary)  and 
circumferential  (secondary)  velocity  at  11  stations  and  also  the 
corresponding  profile  of  the  r.m.s.  turbulent  velocity  profile. 
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The  line  along  which  the  data  are  presented  is  the  radial  line 
bisecting  the  semi-circular  cross  section.  While  this  turned  out 
not  to  be  a  particularly  interesting  line  so  far  as  the  primary 
velocity  was  concerned,  the  secondary  profile  shows  an  interesting 
reversal  near  the  axis  at  135°.  The  breakdown  of  the  classical 
single  cell  that  this  reversal  implies  is  clearly  brought  out  in 
the  computed  secondary  flow  vectors  and  the  corresponding 
streamwise  velocity  contours,  Figure  9. 

While  the  level  of  agreement  obtained  with  the  k~e  model  in 
circular-bend  flows  was  satisfactory  for  many  purposes,  there  was 
nevertheless  room  for  improvement.  Accordingly,  much  of  the  final 
year  has  been  devoted  to  implementing  a  higher  level  turbulence 
model.  While  retaining  the  mixing-length  model  in  the  vicinity  of 
the  wall,  in  place  of  the  k~e  EVM  we  have  incorporated  an  ASM 
version  of  the  Glbson-Launder  [7]  model  for  the  Reynolds  stresses 
(with  minor  adaptations  to  the  model  of  pressure  reflection  from 
walls  to  accommodate  the  toroidal  geometry)  and  for  the  heat 
fluxes  the  generalized  gradient  diffusion  hypothesis: 

—5-  k  -  30 

-  u,0  *  cn  —  u,  u.  -5— 
i  0  e  k  1  dx. 
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An  explanation  of  the  nomenclature  and  further  details  are  given 
in  ref.  [8]  which  appears  as  Appendix  5.  The  more  elaborate  model 
for  the  turbulent  stresses  brought  in  its  train  many  difficulties 
in  securing  convergent  behaviour.  Very  briefly,  the  practices 
that  eventually  proved  most  successful  were: 

(i)  the  evaluation  of  the  Reynolds  stresses  on  a  mesh 
staggered  with  respect  to  the  velocity  nodes  as  shown 
in  Figure  10; 

(ii)  the  explicit  evaluation  of  first  the  normal  stresses 
then  subsequently  the  shear  stresses  over  any  plane 
by  essentially  a  Jacobi  point-substitution  approach. 
(Attempts  at  direct  simultaneous  solution  of  the 
stress  matrix  proved  disastrous); 
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(ill)  the  inclusion  of  parts  of  the  stress  terms  in  the 
momentum  equations  as  effective  viscosities  in  order 
to  promote  diagonal  dominance  of  the  difference 
equations  and  reduce  the  magnitude  of  the  source 
terms : 

(iv)  the  treatment  in  the  turbulence  model  of  negative 
elements  of  the  turbulent  energy  production  rate  as 
dissipation  terms.  (The  somewhat  arbitrary  nature  of 
this  adjustment  is  mitigated  by  the  fact  that  in  the 
final  converged  solution  this  control  is  actually 
invoked  at  only  a  very  small  number  of  nodes.  During 
progress  to  a  converged  field,  however,  the  ASM 
scheme  -  unlike  the  eddy-viscosity  model  employed  in 
our  earlier  work  -  can  allow  large  negative 
generation  to  occur  and  entirely  wreck  what  had 
seemed  to  be  a  smoothly  converging  solution) . 

Appendix  5  gives  a  more  complete  account  of  the  above  practices. 

The  code  with  the  new  model  of  turbulence  incorporated  has 
been  applied  to  re-calculate  the  flow  and  heat  transfer  pattern  in 
the  180°  bend  considered  in  the  experimental  part  of  the  UCB-UMIST 
research.  Extensive  comparisons  are  drawn  in  ref.  [8]  which 
suggest  some  Improvement  in  the  flow  field  predictions  results. 
Figure  11  shows  the  computed  local  Nusselt  number  distribution  at 
several  stations  around  the  U-bend  compared  with  the  experimental 
results  discussed  in  Section  2.  We  see  that  generally  the 
augmentation  of  Nu  on  the  outside  of  the  bend  and  its  damping  on 
the  inside  is  captured  better  by  the  ASM  model  than  with  the  k  e 
EVM.  However,  at  90°  and  135°  the  predicted  level  of  Nusselt 
number  on  the  outside  of  the  bend  is  still  moderately  below  the 
measured  levels.  Recently  repeated  experiments  at  the  same 
Reynolds  number  suggest  that  levels  of  Nu  near  the  outside  of  the 
bend  measured  in  the  first  test  (which  are  those  shown  in  Figure 
11)  may  be  too  high  by  6-8S.  Moreover,  the  neglect  of  density 
changes  in  the  computations  would,  it  is  estimated,  lead  to  an 
underestimate  of  Nu  by  some  3t.  Thus,  about  one  half  the  dis¬ 
crepancy  between  the  ASM  predictions  and  the  data  can  be 
attributed  to  factors  other  than  shortcomings  in  the  turbulence 
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model . 


5 .  Achievements  and  Overall  Conclusions 


The  research  programme  now  concluded  had  set  out  to  establish 
how  well  current  three-dimensional  flow  solvers  could  predict  heat 
transfer  behaviour  in  a  complex  internal  flow  typical  of 
convective  heat-exchanger  situations.  In  order  to  make  the 
assessment,  an  extensive  data  base  had  to  be  acquired,  for 
measurements  of  the  required  level  of  detail,  including  flow  and 
heat  transfer  data  for  the  same  configurations,  simply  did  not 
exist . 


So  far  as  that  data  base  is  concerned,  we  believe  that,  in 
conjunction  with  the  flow  field  and  turbulence  data  gathered  by 
Professor  Humphrey  and  his  team  at  Berkeley,  a  satisfactory  set  of 
experimental  results  has  been  generated.  The  design  of  heat-flux 
meter  developed  in  the  experimental  research  may  also  be  useful  to 
others  embarking  on  constant-wall-temperature  experiments. 

A  more  complicated  picture  is  presented  on  the  computational 
front.  For  the  circular-sectioned  bend,  rather  satisfactory 
agreement  has  been  obtained  especially  with  the  algebraic  stress 
model  of  turbulence.  The  agreement  is  sufficiently  good  that  one 
could  with  confidence  use  computer  modelling  to  optimize  designs 
(excluding  the  1 ow-Reynolds-number  regime  (Re  <  1.5  x  10  )  which 
has  not  been  explored  in  the  present  work). 

For  the  square  bend,  a  quite  different  picture  emerges. 

While,  with  the  adoption  of  PSL  and  the  low-Re  near-wall 

modelling,  the  first  45°  of  development  is  well  accounted  for  with 

the  k^e  model,  there  remains  a  serious  disagreement  with  the 

o 

measured  velocity  profiles  at  90  and,  to  a  diminishing  extent, 
stations  downstream  therefrom.  Our  computations  have  at  least 
shown  that  numerical  errors  are  not  responsible  for  the  difference 
and  that  the  predicted  flow  pattern  is  highly  sensitive  to  the 
form  of  the  dissipation  rate  equation  used. 


Although  the  equations  look  less  complicated  for  the  square- 


Although  the  equations  look  less  complicated  for  the  square- 
than  the  round-sectioned  bend,  the  former  is  undoubtedly  a  much 
more  complicated  flow  both  physically  and  numerically.  Work  is 
continuing  at  UMIST  on  this  flow  sponsored  by  the  NASA  Lewis 
Laboratories  and  we  shortly  hope  to  be  able  to  report 
fully-converged  ASM  results  for  this  case.  Until  such  predictions 
are  available,  it  is  hard  to  say  whether  the  square-sectioned  bend 
is  just  a  difficult  flow  or  whether  it  is  a  pathological  case.  As 
a  tailpiece.  Professor  Humphrey  (personal  communication)  has 
observed  that  the  differences  between  observation  and  experiment 
could  be  explicable  in  terms  of  the  flow  being  unsteady, 
"sloshing"  back  and  forth  about  the  plane  of  geometric  symmetry. 
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ABSTRACT 


An  instrument  for  the  measurement  of  heat  flux  from  a 
surface  with  a  nearly  uniform  temperature  is  described. 
This  instrument  contains  a  thin  film  electrical  resistance 
heater  imbedded  in  a  copper  cone  which  is  thermally  isolated 
from  the  surrounding  walls.  A  differential  thermocouple 
between  the  copper  cone  and  the  wall  is  nulled  such 
that  the  electrical  power  becomes  a  direct  measure  of  the 
surface  heat  flux.  This  instrument  is  modular  and  can  be 
flush  mounted  in  an  external  surface  or  mounted  in  the  wall 
of  a  duct.  It  provides  high  accuracy  and  has  been  used  to 
measure  the  heat  transfer  coefficients  inside  a  circular 
duct.  For  these  measurements  an  uncertainty  analysis  is 
presented  which  shows  that  this  instrument  has  an 
uncertainty  of  +/-3.  6  percent  for  a  convective  heat  flux  of 
342  W/sy  m.  The  major  source  of  uncertainty  was  the 
surface  area.  With  improvements  in  the  surface  area 
measurement  the  uncertainty  would  be  ap pr o x ima t el y  +/-2 
percent. 
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INTRODUCTION 

In  studies  of  convective  and  radiative  heat  transfer, 
it  is  necessary  to  determine  both  the  wall  temperature  and 
the  heat  flux  from  a  surface.  If  these  measurements  are  to 
be  useful.  the  thermal  boundary  conditions  must  be  known. 
There  are  two  common  thermal  boundary  conditions.  In  the 
first  boundary  condition  the  surface  heat  flux  is  uniform. 
This  is  often  accomplished  by  electrically  heating  the  wall. 
In  order  to  determine  the  heat  transfer  coefficient  under 
this  condition  the  temperature  distribution  is  measured  with 
local  sensors.  The  other  important  and  commonly  used 
boundary  condition  is  a  uniform  wall  temperature.  In  this 
case  the  local  heat  flux  is  measured  with  a  local  heat  flux 
gage 

Emrich  (1),  in  a  recent  review  of  heat  flux  gages, 
suggests  that  local  heat  flux  gage  designs  are  based  on 
three  different  conceptual  methods  of  measuring  heat  flux: 

1  the  heat  flux  may  set  up  a  temperature  gradient-  the 
gradient  method;  2.  the  heat  flux  may  cause  a  temperature 
transient  in  a  ca 1  or imet er-th e  transient  method;  3.  The 
heat  flux  may  be  balanced  against  a  calibrated  heat 
removal-the  balance  method.  There  are  many  different  heat 
flux  gage  designs  based  on  these  conceptual  methods  (2.  3. 
and  4)  Unfortunately .  most  of  these  designs  do  not  provide 
a  uniform  wall  temperature  since  they  are  based  on  methods  1 
and  2  and  require  either  a  temperaf  re  gradient  or 
temperature  transient  to  provide  the  sensing  input.  For 
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example*  the  Garden  or  Asymptotic  (thin  circular  foil)  heat 
flux  meter  uses  the  temperature  difference  between  the 

center  and  edge  of  a  surface  mounted  disc  as  the  sensing 
input  (5  and  6)  which  results  in  a  temperature  distribution 
on  the  surface  (method  1-gradient).  This  gradient  in  the 
surface  temperature  may  change  (perturb)  the  surface  heat 
flux  from  that  which  would  occur  in  the  absence  of  the 

sensor.  This  has  been  studied  extensively  for  the  case  of 
convective  heat  transfer  (6-9).  Another  type  of  heat  flux 

gage  using  method  1*  the  wafer-type  meter  (10).  senses  the 

temperature  gradient  through  a  layer  (wafer)  on  a  surface  as 
a  measure  of  the  heat  flux.  This  changes  the  local  surface 
temperature  from  what  it  would  be  in  the  absence  of  the 
sensor  and  thereby  may  change  the  local  heat  flux. 
Thermopile  sensors  (using  differential  thermocouples)  may  be 
fabricated  with  different  geometries  and  be  similar  1 
either  the  asymptotic  or  wafer-type  meters. 

A  heat  flux  gage  which  uses  method  2  (transient)  is  the 
slug-type  heat  flux  meter.  The  slug-type  has  a  thermally 
isolated  mass  and  measures  its  temperature  response  to 
determine  the  heat  flux.  This  changing  surface  temperature 
boundary  condition  may  also  influence  the  local  heat  flux 
(7).  Another  example  of  the  transient  method  (method  2)  is 
the  thin-film  heat  flux  sensor  which  uses  a  surface 
temperature  transient  to  infer  the  heat  flux  (11).  These 
latter  types  of  sensors  are  often  used  in  shock  tunnels  with 
high  temperatures  and  short  flow  periods  Under  these 
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specialized  conditions  the  small  initial  temperature  rise  of 
the  surface  mag  not  significantly  perturb  the  local  heat 
flux. 


As  noted,  each  of  the  above  instruments  perturbs  the 
surface  temperature  while  measuring  the  surface  heat  flux 
which  may  in  turn  perturb  the  surface  heat  flux.  In 
convective  heat  transfer.  this  perturbation  may  persist 
downstream  of  the  sensor  (e.g.,  during  the  development  of  a 
thermal  boundary  layer). 

The  instrument  described  in  reference  12  (which  is 
similar  in  concept  to  that  used  in  reference  13)  is  based  on 
method  3  (the  balance  method)  and  is  designed  to  provide  a 
nearly  uniform  surface  temperature.  In  this  heat  flux  gage 
a  thermally  isolated  but  highly  conducting  mass  (e.  g.  .  a 
copper  cone)  is  mounted  in  a  wall  and  electrically  heated  to 
the  surrounding  surface  temperature.  The  differential 
temperature  between  the  mass  and  the  wall  is  measured  and 
the  power  is  adjusted  until  this  difference  is  reduced  to 
zero  (i.e.  »  nulled)  and  the  null  electrical  power  is  then 
used  to  infer  the  heat  flux. 

In  this  paper,  an  instrument  which  operates  on  the  same 
principle  as  that  of  reference  12,  but  with  several 
important  improvements,  is  described.  This  new  instrument 
is  modular  in  design,  making  it  far  easier  to  fabricate.  It 
uses  a  thin-film  resistor  heater  not  previously  used  in  this 
type  of  sensor  and  unlike  the  earlier  instrument  designs 
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based  on  this  methodi  it  can  also  be  installed  and  removed 
from  a  wall  so  as  to  measure  the  heat  flux  on  an  exterior  or 
an  interior  surface  (i.e.  inside  a  duct). 

DESIGN  OF  THE  HEAT  FLUX  SENSOR 

A  sectional  diagram  of  the  heat  flux  sensor  is  shown  in 
figure  1  and  a  photograph  of  the  sensor  and  its  components 
is  shown  in  figure  2.  The  sensor  is  designed  to  fit  into 
the  wall  of  a  test  rig  as  shown  in  this  figure.  The  face  of 
the  sensor  (top  of  figure  1)  can  be  inside  a  duct  or  on  an 
outside  surface.  The  main  element  of  the  sensor  consists  of 
a  small  copper  cone  with  a  cylindrical  extension  This 
copper  cone  has  an  imbedded  differential  thermocouple  to 
measure  its  temperature  relative  to  the  wall  of  the  test  rig 
and  it  has  a  heating  resistor  to  heat  the  cone  and  balance 
the  surface  heat  flux.  The  copper  cone  is  cemented  into  an 
epoxy  cylinder  which  has  been  drilled  and  machined  to  fit 
the  cone.  The  epoxy  cylinder  fits  into  the  wall  of  the  test 
apparatus  or  surface  of  interest.  This  epoxy  cylinder  has 
an  outside  diameter  of  12  mm  and  has  a  small  groove  machined 
around  the  outside  circumference  with  a  small  access  hole 
for  the  differential  thermocouple  leads  (see  figure  1). 

The  differential  thermocouple  nsists  of  EP-EN-EP 
(Chromel-Constantan-Chromel )  Teflon-covered  wires  (AWG  No 
40  0  076  mm/0.003  in.  diameter)  which  provide  an  output 

of  approx imately  59  microvolts  per  degree  C  The  EP  wire 
enters  from  the  back  of  the  instrument  and  is  then  wound 
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around  the  conical  surface  before  it  is  imbedded  in  the  side 
of  the  conical  surface  where  it  makes  a  junction  with  the  EN 
wire.  The  EN  wire  continues  the  circuit  and  is  also  wound 
around  the  cone  before  it  goes  through  the  small  hole  in  the 
epoxy  cylinder  and  then  around  to  the  opposite  side  of  the 
epoxy  cylinder  where  a  junction  with  the  final  EP  wire  is 
made.  This  final  EP  wire  continues  the  circuit  around 
inside  the  groove  and  enters  the  same  small  side  hole  and 
exits  through  the  back  side  of  the  sensor.  This  arrangement 
considerably  reduces  conduction  between  the  cone  and  outside 
surface  along  the  wires.  It  also  ensures  good  thermal 
contact  between  the  outside  junction  and  the  wall  since 
there  is  a  large  contact  area.  Silver  loaded  epoxy  is 
loaded  into  the  outer  groove  to  further  ensure  that  the 
outside  junction  region  of  the  differential  thermocouple  is 
at  the  wall  temperature.  The  differential  thermocouple 
leads  (EP.EP)  go  to  a  data  logger  where  they  are  connected 
directly  to  an  isothermal  connecting  block.  This  is  shown 
in  the  wiring  diagram  of  figure  3.  In  order  to  protect 
these  delicate  wires  from  mechanical  damage  and  electrical 
interference)  a  copper  screening  braid  is  passed  over  the 
wires  and  anchored  securely  to  the  sensor  cylinder  with  a 
silver  steel  pin.  The  data  logger  that  is  used  to  monitor 
the  signal  from  the  differential  thermocouple  has  a 


resolution  of  10  microvolt. 
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The  heater  which  is  inserted  into  the  copper  cone 
(cylindrical  extension)  is  a  high  stability  metal  film 
heating  resistor  with  a  resistance  of  4.7  Kohms  and  an 

outside  diameter  of  24  mm.  The  active  resistance  part  of 

of  a  ceramic  rod . 

this  resistor  is  a  thin  film  on  the  surface  J  This  has  good 
thermal  contact  with  the  copper  cone  since  the  resistor  is 
cemented  into  the  interior  of  the  cylindrical  extension 
using  silver  loaded  epoxy.  The  electrical  connections  to 
this  resistor  are  made  with  Teflon  coated  constantan  wire  (AUG 
No.  40:  0.076  mm/0.003  inch  diameter)  and  the  ends  are 
coated  with  epoxy  to  avoid  electrical  contact  with  the 
copper  cone.  Four  wires  (two  on  each  end)  are  used  to 
eliminate  lead  wire  resistance  effect  (i.e..  two  wires  carry 
the  current  and  two  wires  are  used  for  the  voltage 
measurement).  In  order  to  minimize  conduction  losses  along 
these  wires-  they  are  routed  so  that  they  first  emerge  from 
the  sensor  body  directly  under  a  copper  spacer  washer  and  in 
good  thermal  contact  with  it  before  being  fed  through  a 
suitable  hole  to  emerge  from  the  rear  of  the  sensor.  These 
fine  wires  are  then  connected  to  four  terminal  posts  secured 
on  the  back  of  the  sensor  where  larger  lead  wires  (0.  56 
mm/0.022  in.  )  are  taken  to  the  power  supply  and  measurement 
instrumentat i on  The  wiring  diagram  for  the  sensor  given  in 
figure  3  shows  the  power  supply  and  measurement  arrangement 
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Following  the  approach  adopted  by  Kraabel  (12)  we  begin 
with  an  energy  balance  equation  for  the  copper  cone  of  the 
heat  flux  sensor. 

Heat  Transfer  Rate  =  Heat  Transfer  Rate  +  Heat  Conduction 
into  the  Cone  from  surface  from  the  Cone 

(a)  (b)  (c) 

ps  =  is  +  g.  01 


The  three  terms  in  the  above  balance  equation  can  be 
expressed  individually  as 


T  7. 

fs  ”  uj  ~ 

fc  + 

(a) 

fcp-s 

(  £rh+  rzai  > 

(b) 

P> 

(c) 

CeCT5'Tw) 

CO 

The  sensor  electrical  input 

p  ower 

is  put  in  terms 

of 

the  sensor  voltage>  shunt  voltage  and  resistance!  and  a  lead 
wire  resistance.  The  lead  wire  resistance  is  only  that 
portion  of  the  lead  wire  in  which  the  ohmic  heating  is 
transferred  by  conduction  to  the  sensor  cone.  It  is 
negligible  in  this  case  since  it  is  much  smaller  than  the 
heater  resistance  of  47  Kohms.  It  is  included  here  for  the 


uncertainty  analysis 
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The  parameter  Ce  is  the  overall  thermal  conductance 
between  the  copper  cone  of  the  sensor  and  the  surrounding 
wall.  Its  value  can  be  determined  experimentally  by 
insulating  the  surface  face  of  the  sensor  and  applying  power 
to  the  heater  resistor.  This  has  been  done  and  the  results 
are  shown  in  figure  4.  The  linearity  of  these  results 
suggests  that  under  these  conditions  the  losses  from  the 
copper  cone  are  primarily  conduction  and  that  the  use  of  Ce 
is  an  appropriate  model.  The  value  of  Ce  can  be  found  from 
the  slope  of  this  curve  and  is  25.7  mW/K. 

Equations  1-4  can  be  combined  to  obtain  an  expression 
for  the  heat  flux  in  terms  of  the  various  measurands. 


This  expression  can  be  used  to  determine  the  individual 
contributions  of  each  measurand  end  the  total  uncertainty  in 
the  measured  heat  flux.  Following  standard  procedure  (see 
reference  2)  the  uncertainty  due  to  each  individual 
measurand  is  given  by 

•»  cJ%. " 

=  ^  (4) 


where  xi  is 
uncertainty 
then  given  by 


the  particular  measurand  and  i  is  the 
in  that  measurand.  The  total  uncertainty  is 


II 
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The  estimated  uncertainties  in  each  measurand.  their 
individual  contribution,  and  the  total  uncertainty  is  shown 
in  Table  I  for  one  set  of  measurements  in  the  application 
described  below  These  estimates  are  based  on  20: 1  odds. 
For  a  measurement  of  342  W/  sq.  m.  the  uncertainty  was 
found  to  be  +/-3.  6  percent.  The  major  source  of  uncertainty 
was  the  surface  area  With  improvements  in  the  surface  area 
measurement  the  uncertainty  would  be  approximately  +  /-2 
percent 

APPLICATION  AND  MEASUREMENT 

Although  this  sensor  can  be  applied  to  many  interesting 
convective  heat  transfer  measurement s,  it  was  developed 
specifically  to  study  convective  heat  transfer  in  pipes  with 
bends  (i.  e.  ,  a  curved  circular  duct)  In  this  application, 
the  face  of  the  sensor  is  machined  with  curvature  in  two 
planes,  which  is  not  possible  with  most  other  heat  flux 
sensors 

In  order  to  test  the  performance  of  the  sensor,  it  was 
first  used  to  measure  heat  flux  rates  in  a  straight  pipe 
with  fully  developed  turbulent  flow  and  with  3  uniform  wall 
temperature.  This  case  was  selected  because  it  has  been 
extensively  studied  and  established  correlations  exist  for 
the  heat  transfer  under  these  conditions. 
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The  test  section  consisted  of  an  externally  heated 
thick  uialled  aluminum  pipe  with  an  internal  diameter  of  7  6 
cm  <3.0  in.  ).  Five  sensors  were  installed  in  the  pipe  at 
equally  spaced  positions  around  the  top  half  circumference 
at  a  position  approx imately  42  diameters  downstream  in  the 
heated  section  The  velocity  and  temperature  profiles  at 
this  position  were  determined  with  a  traverse  probe.  This 
provided  a  measurement  of  the  local  bulk  temperature.  The 
flow  rate  was  determined  using  an  orifice  plate  (constructed 
according  to  British  Standard  1042)  located  far  downstream 
of  the  test  position.  To  ensure  fully  developed  flow 
conditionsi  there  was  an  additional  36  diameters  of  unheated 
pipe  upstream  of  the  heated  section. 


After  the  axisymmetry  of  the  velocity  and  temperature 
fields  in  the  pipe  are  established,  a  series  of  three  tests 
were  made  at  three  different  Reynolds  numbers.  The  heat 
flux  and  corresponding  Nusselt  number  was  then  calculated 
for  each  sensor  as  shown  in  Table  II.  The  Nusselt  number  is 
given  by 


w 


<4 


_ 2. 


Note  that  the  agreement  between  the  different  individual 
sensors  is  extremely  close  for  a  given  Reynolds  number. 
These  results  are  plotted  in  figure  6  and  compared  to  the 
correlation  recommended  by  Petukhov  (seeKays  and  Crawford, 
reference  14)  in  the  form 

.  _  fte  Pr-  ^  (^) 
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In  order  to  determine  if  the  variation  in  measurements  at 
different  circumferential  positions  mas  caused  by  variations 
in  the  individual  heat  flux  sensors.  their  positions  mere 
interchanged  several  times.  It  was  found  that  readings  by 
different  sensors  at  the  same  positions  and  test  conditions 
varied  by  less  than  3‘/.  which  is  within  the  estimated 
uncertainty.  It  was  concluded  that  most  of  the  variation  in 
the  readings  was  due  to  variations  in  the  actual  local  heat 
flux  caused  by  slight  variations  in  the  local  wall 
temperature. 

The  time  response  of  the  sensor  was  investigated  for 
this  application  by  applying  a  step  change  in  the  heater 
voltage  during  a  measurement.  The  resulting  transient 
response  of  the  differential  thermocouple  voltage  can  then 
be  used  to  determine  the  time  constant  for  the  sensor  under 
these  conditions.  The  results  of  a  particular  measurement 
are  shown  in  figure  7.  The  time  constant  under  these 
conditions  was  found  to  be  43  seconds.  This  relatively  long 
time  constant  makes  it  easy  to  null  the  instrument  and  is 
not  limiting  on  the  application  of  the  sensor  in  most  heat 
transfer  investigations  since  steady  conditions  for  most 
experiments  take  much  longer  times. 
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CONCLUSIONS 

A  new  design  for  a  balance  type  heat  flux  sensor  has 
been  described  The  primary  advantage  of  this  design  over 
earlier  sensors  is  its  modular  characteristic  and  ability  to 
be  inserted  into  a  wall  to  study  either  exterior  or  interior 
heat  transfer.  It  also  uses  a  thin  film  resistor  heating 
element  not  used  in  earlier  designs.  Unlike  earlier  sensors 
of  this  type,  it  can  and  has  been  fabricated  in  quantities. 
Tests  with  different  individual  sensors  under  the  same  test 
conditions  show  that  the  sensors  produce  repeatable  results 
within  the  estimated  uncertainty  of  the  measurements. 


ACKNOWLEDGEMENT.  The  authors  gratefully  acknowledge  the 
support  of  the  U.  S.  Office  of  Naval  Research  (Power 
Program)  through  grant  N00014-83-G-0021  under  the  direction 
of  Professor  B  E.  Launder  at  UMIST.  This  c ol 1 ab ora t i ve 
effort  was  also  made  possible  in  part  by  the  USA. 
National  Science  Foundation  grant  MEA-84-41626 

NOTE  Authors  are  listed  in  alphabetical  order 


Page  15 


REFERENCES 

1.  Enrich  R.  J.  "Methods  of  Experimental  Physics.  Vol 
18B  Fluid  Dynamics. "  Chapter  7-Heat  Transfer  Gages  by  W. 
Paul  Thompson,  pp  663-685,  Academic  Press  198 1. 

2.  Holman,  J  P  ,  "Experimental  Methods  for  Engineers, " 

Fourth  Edition,  McGraw  Hill,  1984. 

3.  Doebelin,  E  O  ,  “Measurement  Systems:  Application  and 

Design,"  McGraw  Hill,  1975. 

4.  Eckert,  E  R.  G  .  and  Goldstein,  R.  J.  ,  "Measurements 

in  Heat  Transfer,"  Second  Edition,  McGraw  Hill,  1976. 

5.  Hornbaker,  D  R.  ,  and  Rail,  D.  L.  ,  "Thermal 

Perturbations  Caused  by  Heat  Flux  Transducers  and  Their 
Effect  on  the  Accuracy  of  Heating  Rate  Measurements. "  ISA 
Transac t i ons ,  p  100,  April  1963. 

6.  Garden,  R  ,  "An  Instrument  for  the  Measurement  of 
Intense  Thermal  Radiation,"  Review  of  Scientific 
Instruments,  Vol.  24,  No.  5,  1953,  pp.  366-370. 

7.  Gardon,  R.  ,  "A  Transducer  for  the  Measurement  of  Heat 

Flow  Rate,"  ASME  Tr an  sac t i on s ,  1960,  pp.  396-398. 

8  Hager,  N.  E  ,  Jr  ,  "Thin  Foil  Heat  Meter, "  Review  of 

Scientific  Instruments,  November  1965. 

9  Malone,  E.  W  ,  "Design  and  Calibration  of  Thin-Foil 

Heat  Flux  Sensors, "  ISA  Transac ti ons,  Vol.  7,  No  3,  pp 
175-180 

10  Streigl.  S  A.  ,  and  Diller,  T.  E.  ,  "The  effect  of 

Entrainment  Temperature  on  Jet  Impingement  Heat  Transfer, " 
ASME  Journal  of  Heat  Transfer,  Vol  106,  1984,  pp.  27-33. 

11.  Cook,  W.  J.  ,  "Determination  of  Heat  Transfer  Rates 

from  Transient  Surface  Temperature  Measurement, "  AIAA 
Journal.  Vol.  8  .  No.  7,  pp.  1366-1368. 

12.  Kraabel,  J  S. ,  Baughn,  J.  W  .  and  McKillop,  A.  A. 

"An  Instrument  for  the  Measurement  of  Heat  Flux  from  a 
Surface  with  Uniform  Temperature."  ASME  Journal  of  Heat 
Transfer,  Vol  102.  PP576-578,  August  1980 

13.  Achenbach.  E  .  "Total  and  Local  Heat  Transfer  from  a 
Smooth  Circular  Cylinder  in  Cross-Flow  at  High  Reynolds 
Number, "  Int  J  Heat  Mass  Transfer. ,  Vol.  18,  pp 
1387-1396 

14.  Kays,  W  M  ,  and  Crawford,  M  E.  ,  Convective  Heat  and 
Mass  Transfer,  McGraw  Hill,  1980,  pg  243 


Page  16 


NOMENCLATURE 


-RH 


-RS 


RS 


!rh 

K 

Nu 

Pr 

Ps 

9c 


RS 

Re 

Rlw 

TB 


sensor  face  area  (n  ) 

sensor  thermal  conductance  (W/*C) 

friction  coefficient 

heater  resistor  voltage  (V) 

shunt  resistor  voltage  (V) 

differential  thermocouple  voltage  ( uV ) 

shunt  resistor  current  (ft) 

heater  resistor  current  (ft) 

fluid  thermal  conductivity  (W/m*C) 

Nusselt  number 

Prandtl  number 

sensor  power  (W) 

conduction  neat  rate  (W) 

sensor  surf  ace  heat  rate  (W) 

2 

sensor  heat  flux  (W/m  ) 
shunt  resistance  (fi) 
heater  resistance  (fi) 

Reynolds  number 
lead  wire  resistance  (fi) 
fluid  bulk  temperature  (*C) 
sensor  temperature  (*C) 
wall  temperature  (’C) 
independent  variable 


TABLE  I 
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xi 

Value  of 

xi 

4Xi 

M 

Sqs  «x.  (w / m2 ) 
axi 

(%) 

ers 

1.3V 

0.005V 

1.33 

0.37 

erh 

13.00V 

0.05V 

1.33 

0.37 

«H 

4700  n 

l.o  n 

.73 

0.21 

R£w 

11.9  n 

1.0  0 

.073 

.02 

As 

105.irm2 

3.5mn? 

11.42 

3.33 

VTs 

0*K 

.017*K 

4.15 

1.21 

TABLE  II 

NUSSELT  NUMBER  MEASUREMENTS  WITH 

HEAT 

FLUX  SENSORS 

RE 

NOi 

nu3 

* 

nu4 

nu5 

33,700 

75.33 

73.31 

74.43 

70.60 

71.45 

58,600 

121.6 

122.1 

125.8 

116.8 

118.2 

87,500 

165.5 

162.5 

_ 

159.6 

162.1 

Diagram  of  the  Heat  Flux  Sensor 


Figure  3  -  Instrumentation 
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Figure  4  -  Effect  of  Heater  Power  on  the  Differential  Thermocouple 
Voltage  During  an  Insulated  Test 
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Figure  6  -  Nusselt  Number  Results  In  a  Heated  Duct  Experiment 
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Figure  7  -  Response  of  the  Differential  Thermocouple 
Voltage  to  a  Step  Heater  Voltage  Change 
During  a  Heated  Duct  Experiment 
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Abstract 


The  paper  reports  extensive  convective  heat  transfer  data  for  turbulent 
flow  of  air  around  a  U-bend  with  a  ratio  of  bend  radius:pipe  diameter 
of  3.375  :  1.  Experiments  cover  Reynolds  numbers  from  2  x  10**  -  1.1  x 
10^.  Measurements  of  local  heat  transfer  coefficient  are  made  at  six 
stations  and  at  five  circumferential  positions  at  each  station.  At 
Re  =  6  x  10^  a  detailed  mapping  of  the  temperature  field  within  the  air 
is  made  at  the  same  stations.  The  experiment  duplicates  the  flow 
configuration  for  which  Azzola  and  Humphrey  [3]  have  recently  reported 
laser-Doppler  measurements  of  the  mean  and  turbulent  velocity  field. 

The  measurements  show  a  strong  augmentation  on  the  outside  of  the  bend 
and  relatively  low  levels  on  the  inside  associated  with  the  combined 
effects  of  secondary  flow  and  the  a mpl i f i ca t i on / supp r e s s ion  of 
turbulent  mixing  by  streamline  curvature.  The  peak  level  of  Nu  occurs 
half  way  around  the  bend  at  which  position  the  heat  transfer 
coefficient  on  the  outside  is  about  three  times  that  on  the  inside. 
Another  feature  of  interest  is  that  a  strongly  non-uniform  Nu  persists 
six  diameters  downstream  of  the  bend  even  though  secondary  flow  and 
streamline  curvature  are  negligible  there.  At  entry  to  the  bend  there 
are  signs  of  partial  1  a m i n a r i z a t i o n  on  the  inside  of  the  bend,  an 
effect  that  is  more  pronounced  at  lower  Reynolds  numbers. 
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1  .  Introduction 


Turbulent  flow  in  a  pipe  bend  is  fundamental  to  many 
engineering  components  in  the  thermal  power  area,  especially  those 
associated  with  convective  heat  transport.  As  two  examples  may  be 
mentioned  the  liquid  side  of  a  shell-and-tube  heat  exchanger  and 
the  interior  cooling  passages  of  turbine  blades.  The  flow 
geometry  is  moreover  sufficiently  basic  that  -  even  for  laminar 
flow  -  it  has  been  the  subject  of  numerous  experimental, 
analytical  and  numerical  studies. 

In  the  turbulent  flow  regime  the  extensive  pitot  tube 

explorations  by  Rowe  [1]  of  the  approach  to  its  fully  developed 

state  of  flow  in  a  gently  curved  U-bend  threw  considerable  light 

on  the  markedly  different  flow  structure  at  different  positions  in 

the  bend.  More  recently  two  LDA  explorations  have  helped  provide 

more  extensive  mappings  of  the  flow  structure.  Enayet  et  al  [2] 

measured  the  distribution  of  streamwise  mean  and  rms  fluctuating 

velocity  along  several  lines  traversing  the  duct  cross  section  at 

five  stations  within  and  downstream  of  a  90°  bend  with  mean  bend 

radius  2.8  times  the  pipe  diameter.  The  mean  velocity  profiles 

showed  clearly  the  piling  up  of  low  momentum  fluid  on  the  inside 

of  the  bend  by  the  secondary  flow.  This  low  momentum  fluid  was 

further  retarded  by  the  adverse  pressure  gradient  the  flow  on  the 

inside  encountered  as  it  passed  from  the  curved  bend  to  the 

straight  exit  tangent.  The  study  by  Azzola  and  Humphrey  [31 

obtained  measurements  of  the  mean  and  rms  values  of  both  the 

streamwise  and  the  circumferential  (secondary)  components  along 

the  diametral  line  at  right  angles  to  the  duct's  plane  of 

symmetry.  The  data  were  obtained  for  two  turbulent  Reynolds 

U  S 

numbers  (Re  .  =  5.79  x  10  and  1.1  x  10)  and  at  a  total  of  12 

d 

stations  including  five,  at  one-diameter  intervals,  immediately 
downstream  of  the  bend.  Their  measurements  of  the  secondary  flow, 
summarized  in  Figure  9  and  discussed  more  fully  in  Section  3.  show 
in  quantitative  detail  a  number  of  the  qualitative  inferences 
drawn  by  Rowe  [1],  including  the  reversal  of  the  secondary  flow  in  ; 

the  centre  of  the  pipe  by  135°  into  the  bend. 

t 

No  correspondingly  detailed  sets  of  data  for  heat  transfer  in  i 

c i r cu 1  a r -se c t i on ed  pipe  bends  are  available.  Local  heat-transfer  , 

1  \ 
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coei i 1 cients  for  turbulent,  nominally  fully-developed  flow  have 
been  reported  by  Seban  and  McLaughlin  [4].  Their  experiments, 
carried  out  over  a  range  of  Reynolds  number,  showed  mean  levels  of 
h e a t - t r a n s f e r  coefficient  some  505  above  those  for  a  straight  tube 
and  of  the  order  of  2:1,  the  levels  diminishing  mono ton  i  c a  1 1  y  from 
the  outside  to  the  inside  of  the  bend.  Qualitatively  this 
behaviour  is  to  be  expected:  the  secondary  flow  continually 
removes  hot  fluid*  from  the  outside  of  the  bend  and  transports  it 
to  the  inside.  There  is  a  further  contributor  to  the  circum¬ 
ferential  non-uniformity.  Streamline  curvature  will  tend  either 
to  dampen  or  augment  turbulent  mixing  depending  on  whether  the 
mean  angular  momentum  increases  or  decreases  with  radius  (here 
referenced  to  the  origin  of  the  streamline,  not  the  pipe  centre), 
Bradshaw  [5].  This  effect  reduces  thermal  mixing  on  the  inside  of 
the  bend  and  increases  it  on  the  outside,  i.e.  its  effect  on  the 
heat  transfer  coefficients  in  the  same  sense  as  that  of  the 
secondary  flow. 

The  aim  of  the  experimental  work  presented  in  this  paper  has 
been  to  provide  detailed  information  about  the  convective  heat 
transfer  coefficient  and  thermal  field  in  a  developing  flow  around 
a  bend,  a  more  complicated  and  more  important  situation  from  a 
practical  standpoint  than  fully-developed  flow.  Because  of  the 
desirability  of  having  flow-field  data  as  well  (since  convection 
by  the  secondary  flow  and  modifications  to  the  fluctuating 
velocity  are  the  agencies  determining  the  level  of  heat  transfer 
coefficient),  the  experiment  has  been  designed  to  replicate  the 
geometry  and  boundary  conditions  of  the  180°-bend  experiments  of 
Azzola  and  Humphrey  [ 3 ] .  Together,  those  experiments  and  the 
present  ones  provide  the  corresponding  mapping  for  the  circular 
tube  that  the  recent  experiments  by  Chang,  Humphrey  and  Modavi  [6] 
and  Johnson  and  Launder  [7]  have  reported  for  the  square-sectioned 
U-bend  . 
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assuming  a  heated  tube  wall 
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2 .  Apparatus.  Instrumentation  and  Experimental  Procedure 

The  test  section,  shown  in  Figure  1,  consisted  of  a  76mm 
diameter  c i r c u 1  a r -s e c t i o n e d  U-bend  with  inlet  and  exit  tangents 
each  50  diameters  long,  25  of  which  were  instrumented  for  heating 
(see  Figure  1).  A  short  perspex  spacer  was  interposed  between  the 
flanges  separating  the  heated  and  unheated  halves  of  the  inlet/ 
exit  tangents  to  keep  heat  conduction  to  the  unheated  sections  to 
unimportant  levels.  The  apparatus  was  designed  to  provide  a 
nearly  uniform  wall-temperature  boundary  condition,  the  tangent 
tubing  being  made  from  sections  of  thick-walled  drawn  aluminium 
tubing  while  the  180°  bend  itself  was  machined  in  two  symmetrical 
halves  from  solid  aluminium  blocks.  The  mean  radius  of  the  bend 
was  6.75  times  the  pipe  radius.  The  sections  of  the  inlet  and 
outlet  tangents,  each  0.95m  in  length,  were  flanged  and  spigoted 
at  their  ends  to  facilitate  their  connection  together  without  mis¬ 
match  at  the  pipe  wall;  for  the  same  reason  the  upper  and  lower 
halves  of  the  180°  bend  were  also  spigoted. 

The  flow  entered  the  test  section  via  a  16:1  area 
contraction.  The  flow  in  the  inlet  tangent  was  rather  sensitive 
to  the  flaring  of  the  bellmouth  and  some  of  the  earlier  tests  were 
taken  without  the  contraction  as  this  led  to  a  more  symmetrical 
velocity  distribution.  Later  tests,  with  a  re-designed  bellmouth 
in  place,  achieved  a  further  improvement  in  symmetry.  Air  was 
drawn  through  the  rig  by  a  10  HP  blower  which  provided  a  maximum 
test  section  Reynolds  number  of  1.1  x  10'’.  Flow  rates  were 
monitored  by  an  orifice  plate  located  in  a  separate  section 
leading  to  the  blower;  the  plate  itself  and  associated  pipework 
conformed  with  BS1042. 

Thermal  measurements  were  made  at  five  stations  in  the  bend 
(15°,  55°,  90°,  125°  and  165°)  and  at  six  diameters  downstream. 
The  data  comprised  local  heat  transfer  rates  measured  at  five 
circumferential  positions  in  the  top  half  of  the  bend.  These  were 
measured  by  heat-flux  gauges  developed  at  UMIST  from  an  earlier 
design  of  Kraabel,  Baughn  and  McfCillop  [8]  and  described  in  detail 
in  Baughn  et  al  [9].  Briefly,  the  copper  cone  providing  the 
heating  surface,  with  a  small  resistor  embedded  at  its  centre,  was 
enveloped  by  a  cylindrical  body  of  low  thermal  conductivity  epoxy 


(Figure  2b).  The  unit  was  flush  mounted  with  the  surface  of  the 
pipe  bend.  Heat  input  to  the  resistor  was  adjusted  until  the  cone 
temperature  was  brought  to  that  of  the  surrounding  pipe  (to  within 
0.015°C).  Accurate  measurement  of  the  current  supply  to  the 
resistor  and  the  voltage  applied  across  it  enabled  the  heat  flux 
through  the  exposed  area  of  the  sensor  to  be  calculated.  Prior  to 
their  use  in  the  present  study,  the  gauges  were  extensively  tested 
in  a  straight  pipe  [93. 

At  one  Reynolds  number  temperature, distributions  within  the 
air  stream  were  made  by  means  of  a  thermocouple  rake  consisting  of 
seven  PTFE  insulated  chromel-alumel  thermocouples  of  0.2mm 
diameter,  Figure  2c.  These  were  held  by  a  tufnol  (low  thermal 
conductivity)  support  bar  curved  to  match  the  radius  of  the  pipe, 
thus  allowing  all  seven  thermocouples  to  record  temperatures  close 
to  the  wall.  The  rake  was  traversed  vertically  covering  the  lower 
half  of  the  bend.  Because  of  the  time-consuming  nature  of  this 
mapping,  other  tests  merely  recorded  -  with  a  single  thermocouple 
-  the  temperature  at  the  bend  centre.  Unlike  the  rake,  the  single 
probe  could  be  inserted  and  removed  through  the  stem  access  hole 
without  disassembling  the  bend. 

The  heat  supply  to  the  bend  itself  and  to  the  inlet  and  exit 
tangents  was  provided  by  self-adhesive  heat  tape  (Clayborn  Labs. 
Inc.)  affixed  to  the  outside  of  the  bend  with  approximately  12mm 
spacing  between  adjacent  tracks.  The  tape  on  the  straight 
sections  was  wound  helically  while  for  the  bend  section  twenty 
single  tracks  running  axially  were  applied  to  the  upper  and  lower 
halves  of  the  bend.  The  electrical  connections  allowed  the 
heating  current  to  twenty  tracks  on  the  outer  half  of  the  bend  (1C 
each  from  the  upper  and  lower  halves)  to  be  controlled  separately 
from  the  twenty  on  the  inside.  A  two-dimensional  finite- 
difference  he  a t -con d uc t  i  on  analysis  indicated  that  at  the  inside 
surface  of  the  tube  temperature  variations  arising  from  the 
discrete  heating  on  the  outside  amounted  to  no  more  than  2*  of  the 
temperature  difference  between  the  wall  and  ambient  values.  There 
was  found  to  be  a  larger  variation  in  wall  temperature  in  the 
streamwise  direction  (  ~-7%)  arising  from  changes  in  both  the  bulk 
temperature  and  the  velocity  field  as  the  air  passes  around  the 
bend.  This  variation  was  of  course  taken  into  account  in 
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evaluating  the  Nusseit  number  data  presented  in  Section  3. 

Before  the  test  programme  began,  the  apparatus  was  checked 

for  leaks  by  sealing  off  the  pipework  and  pressurizing  the  rig. 
Thereafter  velocity  profiles  were  taken  at  two  diameters  upstream 
of  the  bend  to  establish  that  the  flow  was  acceptably  symmetric 
and  to  provide  well-defined  entry  conditions.  Figure  3  shows  the 
mean  and  rms  axial  velocity  profiles  for  both  vertical  and 
horizontal  traverses;  the  horizontal  traverse  already  shows  a 
shift  of  the  airflow  towards  the  inside  of  the  bend  induced  by  the 
upstream  effects  of  the  bend  on  the  pressure  field.  To  verify 
that  the  cold-rig  flow-field  data  (obtained  by  the  Berkeley  group) 
would  be  representative  of  flow  conditions  in  the  heated  duct,  the 
mean  profile  was  measured  again  with  a  pitot  tube  with  the  duct 

both  heated  and  unheated.  No  significant  differences  in  the 

profile  were  noted  with  the  pipe  wall  temperature  20°C  above  that 

of  the  entering  fluid  both  showing  complete  agreement  with  the 
cold-flow  profile  obtained  by  the  hot  wire. 

The  heat  transfer  measurements  were  taken  at  three  nominal 

HU  U 

Reynolds  numbers;  2  x  10  ,  6  x  10  and  11  x  10  .  The  two  higher 
Reynolds  numbers  correspond  closely  with  those  for  which  the 
velocity  field  has  been  mapped  by  Azzola  and  Humphrey  [33.  In  our 
work  the  greatest  attention  has  been  given  to  the  middle  Reynolds 
number.  Data  have  been  obtained  with  the  nominal  wall  temperature 
10°C  and  20°C  above  that  of  the  air  at  inlet.  Iacovides  [10]  has 
provided  a  detailed  analysis  of  experimental  uncertainties.  Based 
on  this  we  estimate  that  the  Nusseit  number  data  are  accurate  to 
within  ±R%  at  the  highest  Reynolds  number  with  20°C  overheat 
(where  uncertainties  are  least)  and  to  within  ±10%  for  a  10°C 
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overheat  at  Re  =  2  x  10  (where  uncertainties  are  greatest). 


3 .  Presentation  and  Discussion  of  Results 


First  we  re-examine  the  mean  velocity  field  data  of  Azzola 
and  Humphrey  [33  summarized  in  Figure  U  to  provide  a  background 
for  presenting  the  thermal  measurements.  The  velocities  were 
taken  along  the  vertical  diametral  line  perpendicular  to  the  plane 
of  geometric  symmetry.  Fully-developed  pipe  flow  at  entry  leads 
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to  a  rapid  development  of  secondary  velocities  which,  by  **5°  into 
the  bend,  reach  a  peak  value  of  about  ?  0  *  of  the  bulk  streanwise 
velocity,  the  maxinun  occurring  very  close  to  the  wall.  The 
secondary  notion  carries  slow  noving  fluid  near  the  pipe  wall  to 
the  inside  of  the  bend;  cor  respond i ngl y  faster  moving  fluid  close 
to  the  symmetry  plane  migrates  towards  the  outside  of  the  bend 
giving  the  classical  single-cell  vortex.  By  90°  into  the  bend, 
the  convection  of  fluid  by  the  secondary  motion  shifts  the  peak  in 
streamwise  velocity  near  the  pipe  wall  with  a  minimum  at  the  pipe 
centre.  Because  of  this  inversion  the  radial  pressure  gradient 
now  acts  to  drive  fluid  near  the  pipe  centre  to  the  inside  of  the 
bend  (135°)  and  this  reverse  stirring  makes  the  streamwise 
velocity  virtually  uniform  along  the  measurement  line  (177°  to 
four  diameters  downstream). 

Figure  5  provides  the  corresponding  temperature  mapping 

obtained  in  the  present  explorations  at  a  nominal  bulk  Reynolds 

4  o 

number  of  6  x  10  and  for  the  tube  wall  heated  =  20  C  above 

ambient.  The  variation  of  (T  -  T)/(T  -  T.)  along  seven  vertical 

w  /  w  1 

lines  is  shown;  the  centre  line  is  the  same  as  that  on  which  the 
velocity  profiles  of  Figure  U  were  obtained.  Here  T w  is  the  wall 
temperature,  T  the  local  air  temperature  and  T  ^  the  air 
temperature  at  inlet  to  the  pipe.  At  15°  into  the  bend  the 
temperature  field  differs  little  from  that  found  in  a  straight 
pipe:  along  each  traverse  line  the  normalized  temperature  rises 

monoton i ca 1 1 y  from  the  wall  to  the  symmetry  plane  while  along  the 
symmetry  plane  itself  the  temperature  is  nearly  symmetric  with 
respect  to  the  pipe  centre. 

A  sharply  different  picture  has  developed  by  55°.  Along  the 
two  traverse  lines  closest  to  the  inside  of  the  bend,  a  deep 
trough  in  normalized  temperature  appears  in  the  vicinity  of  the 
symmetry  plane  associated  with  hot  air  from  the  immediate 
proximity  of  the  bend  wall  being  lifted  away  from  the  surface  and 
carried  towards  the  pipe  centre  by  the  secondary  flow.  Along  the 
third  line  there  is  only  a  weak  dip  in  normalized  temperature  at 
this  station  and  none  at  all  along  the  fourth  (the  diametral 
line):  this  is  because  at  these  locations  the  secondary  motion 

has  brought  fluid  at  only  a  slightly  different  temperature  than 
that  originally  occupying  that  position.  Over  the  outer  half  of 
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pipe  a  reverse  picture  emerges:  here  the  secondary  flow  "imports" 
faster  moving,  cooler  fluid  (i.e.  with  larger  normalized 
temperatures  )  . 

The  temperature  field  at  90°  is  very  much  a  continuation  of 
the  developments  seen  at  55°.  There  is  now  a  marked  dip  in 
temperature  along  the  centre  traverse  line,  a  feature  which 
matches  the  dip  in  streamwise  velocity  noted  above  in  Figure  9, 
The  profiles  at  125°  and  165°,  like  the  corresponding  velocity 
profiles,  now  show  little  variation  along  the  traverse  lines 
though  there  remains  a  significant  rise  in  (normalized) 
temperature  from  the  inside  to  the  outside  of  the  bend.  By  six 
diameters  downstream  of  the  bend  the  re-establishment  of  new 
thermal  boundary  layers  on  the  pipe  wall  is  evident. 

The  overall  picture  that  emerges  from  the  temperature-field 
data  is  one  that  is  fully  compatible  with  the  A z z o 1  a - H u m ph r e y 
velocity  measurements .  Indeed,  because  the  present  measurements 
provide  a  view  of  the  temperature  field  over  the  whole  cross 
section,  they  allow  a  f.uller,  more  quantifiable  impression  than 
could  be  gained  from  the  flow-field  data  alone. 

The  measured  heat-transfer  coefficients,  expressed  in  terms 
of  Nusselt  number,  are  presented  in  Figure  6.  It  is  emphasized 
that  in  the  definition  of  Nu  the  temperature  difference  is  that 
between  the  local  wall  value  and  the  air  temperature  at  the  centre 
of  the  tube.  This  basis  was  adopted  because  the  bulk  mean 
temperature  at  any  station  could  not  be  determined  with  great 
precision  since  the  experiments  were  with  a  (nearly)  uniform  wall 
temperature.  Figure  6  relates  to  precisely  the  test  conditions 
for  which  the  interior  temperature  profiles  were  obtained 
indeed,  the  rake  was  still  in  place,  pushed  down  as  far  as 
possible  while  the  measured  heat  fluxes  were  obtained  around  the 
top  half  of  the  duct.  Perhaps  the  first  thing  to  notice  is  that 
even  by  15°  an  appreciable  circumferential  variation  of  heat  flux 
has  developed.  Much  stronger  circumferential  variations  are  found 
as  one  proceeds  around  the  bend,  a  development  due  partly  to  the 
secondary  flow  (the  'impingement'  of  cool  fluid  on  the  outer  wall 
will  obviously  raise  heat  transfer  rates)  and  partly  to  the 
modifications  to  the  turbulent  mixing,  augmentation  occurring  on 
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the  concave  surface  and  damping  on  the  inside  of  the  bend.  The 
maximum  mean  level  of  Nu  occurs  at  90°.  The  variation  from  the 
inside  to  the  outside  of  the  bend  is  also  largest  there:  extra¬ 
polation  to  0°  suggests  an  out e r -t o-i nne r  variation  of  3:1.  The 
mean  level  and  the  circumferential  variation  of  Nusselt  number 
falls  gradually  over  the  second  half  of  the  bend  due  to  decay  of 
secondary  flow  (see  Figure  4b).  Even  six  diameters  downstream, 
however,  there  is  still  a  2:1  variation  in  Nu  between  (what  had 
been)  the  inside  and  outside  of  the  bend.  The  reason  for  this 
seems  to  be  that  due  to  the  virtual  disappearance  of  the  secondary 
flow  the  hottest  fluid,  which  during  passage  around  the  bend 
accumulated  on  the  inside,  simply  stays  there  (Figure  5f)  until 
eventually  the  (relatively)  slow  process  of  turbulent  diffusion 
restores  axial  symmetry  to  the  temperature  field.  Evidently, 
however,  this  occurs  only  far  downstream  of  the  region  examined  in 
the  present  experiments. 

More  extensive  heat-flux  measurements  in  the  bend  itself  were 
made  in  which  just  the  centreline  temperature  was  measured;  in 
this  way  the  risk  of  affecting  the  surface  heat  flux  pattern  with 
the  thermocouple  rake  was  avoided.  Figure  7  shows  the  distrib¬ 
utions  of  Nusselt  number  for  nearly  the  same  Reynolds  number  as 
Figure  6  for  two  temperature  differences:  20 °C  (as  in  Figure  6) 
and  1 0°C .  In  a  few  positions  the  heat-flux  guage  went  'open 
circuit'  during  the  test,  which  is  why  data  points  are 
occasionally  missing.  The  variation  is  very  similar  to  that  shown 
in  Figure  6,  though  on  the  whole  the  levels  are  lower  by,  on 
average,  5%.  For  the  test  in  which  data  with  the  temperature 
difference  was  reduced  to  10 °C  the  levels  of  Nu  lie  on  average  8* 
below  those  for  a  20°C  temperature  difference.  Two  earlier  tests, 
in  which  just  the  heat  transfer  rates  were  measured,  also  showed 
the  same  effect,  i.e.  the  heat  fluxes  divided  by  the  temperature 
difference  at  inlet  was  8-105  less  at  the  lower  heating  rate. 

The  decrease  in  density  associated  with  heating  could  be 
expected  to  raise  the  Nusselt  number  in  a  curved  pipe  by  two 
mechanisms.  Firstly,  the  secondary  flow  driving  force  is  raised 
since  the  lowest  density  fluid  is  also  the  slowest  moving. 
Secondly,  there  will  be  some  tendency  for  cooler  (denser) 
turbulent  eddies  to  centrifuge  to  large  radii  thus  promoting 


nixing  on  the  outside  of  the  bend  (though  dampening  it  on  the 
inside).  Neither  mechanism  looks  powerful  enough  to  have  a 
significant  effect,  however,  when  such  small  density  differences 
are  involved.  We  note,  moreover,  that  at  a  Reynolds  number  of 

C 

1.1  x  10  ,  for  which  results  are  shown  in  Figure  8,  there  are  no 

significant  differences  between  the  20°C  and  10°C  experiments. 

Now,  at  this  higher  Reynolds  number  precision  is  greater  because 

heat  input  rates  are  greater  and  heat  losses  of  smaller 

significance.  Our  view  is  thus  that  the  apparent  sensitivity  to 

U 

AT,  shown  in  the  data  at  Re  =  5.8  x  10  ,  is  probably  not  a  real 

effect.  The  values  obtained  with  the  larger  temperature 

difference  are  thus  preferred. 

Finally,  Figure  9  shows  the  results  of  consolidating  data  for 

three  Reynolds  numbers  -  the  two  for  which  data  are  shown  above 

and  for  20,000.  In  all  cases  the  measured  Nusselt  number  is 
0  8 

divided  by  Re  *  for  the  test  in  question.  While  the  data  for  the 
three  tests  broadly  collapse  together,  there  is  evidently  a 
diminution  in  the  ratio  of  outerrinner  heat  transfer  coefficient 
as  the  Reynolds  number  is  raised,  due  mainly  to  a  marked  increase 
in  Nusselt  number  on  the  inside  of  the  bend  -  or,  conversely,  a 
sharp  fall  as  Re  is  lowered.  The  indication  seems  to  be  that  the 

diminished  turbulent  mixing  near  the  inner  wall  is  causing  the 

viscous  sublayer  to  become  thicker  at  the  lowest  Reynolds  number 
causing  a  substantial  increase  in  thermal  resistance.  This  is 

qualitatively  similar  to  the  effects  found  in  other  boundary  layer 
flows  subject  to  damping  of  turbulence.  Indeed,  for  flows  in 
coils  it  is  known  (see,  for  example,  Sreenivasan  and  Strykowski 

H 

[11])  that  at  Reynolds  numbers  below  10  the  flow  on  the  inside  of 

the  bend  may  become  laminar  even  though  turbulent  flow  is 

maintained  on  the  outside. 


Concluding  Remarks 

A  detailed  mapping  has  been  provided  of  the  thermal  field  in 
turbulent  flow  around  a  circular-sectioned  U-bend  for  conditions 
corresponding  to  those  for  which  Azzola  and  Humphrey  have  reported 
velocity  field  data.  Together,  the  two  experiments  give  a 
comprehensive  documentation  of  the  processes  of  convective  heat 
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transport  in  a  conplex  but  well-defined  three-diuensional 
turbulent  flow.  It  is  hoped  that  they  will  provide  a  searching 
test  case  for  the  t h r ee-d i men s i ona  1  flow  computation  procedures 
that  are  now  becoming  widely  available  for  industrial  design 
purposes  . 

Among  the  most  striking  results  to  emerge  were: 

•the  large  circumferential  variation  of  heat  transfer 
coefficient  that  developed  by  only  15°  into  the 

i) 

bend:  nearly  a  2:1  ratio  for  fie  =  2  x  10  ; 

•the  attainment  of  maximum  mean  heat  transfer 
coefficient  at  90°  into  the  bend  at  which  station 
the  maximum  c i r c urn f e r e n t i a  1  variation  was  also 

H 

recorded  (nearly  4:1  for  Re  =  2  x  10  ); 

•the  persistence  of  strong  circumferential  non¬ 
uniformities  in  Nusselt  number  six  diameters 
downstream  of  the  bend. 
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ABSTRACT 


laser-Doppl er  measurements  of  the  longitudinal  and  ci rcumferent i al  velocity 
components  are  reported  for  developing  turbulent  flow  in  a  strongly  curved  180° 
pipe  and  its  downstream  tangent.  In  the  bend,  the  mean  longitudinal  velocity 
component  changes  little  after  8  =  90°,  but  the  circumferential  component  never 
achieves  a  fully-developed  state.  Similar  behavior  is  observed  in  the  normal 
stresses,  with  large  levels  of  flow  anisotropy  arising  everywhere  in  the  bend 
and  downstream  tangent.  Between  9  =  90°  and  X/D  =  5,  the  circumferential  velo¬ 
city  profiles  display  reversals  of  the  secondary  flow  which  are  essentially 
independent  of  the  Reynolds  number. 

Predictions  of  the  flow  development  are  presented  based  on  a 
"semi-el  1 iptic"  truncation  of  the  Reynolds  equations  in  the  main  part  of  the 
flow  with  the  standard  k-e  effective  viscosity  model  used  to  approximate  the 
turbulent  stress  field.  In  the  immediate  vicinity  of  the  wall  a  simpler  treat¬ 
ment,  PSl,  is  adopted  that  allows  the  inclusion  of  the  very  fine  mesh  needed  to 
resolve  the  viscous  sublayer  without  excessive  computer  storage.  The  calculated 
behavior  displays  reasonably  good  agreement  with  the  measurements  in  the  bend, 
including  the  secondary  flow  reversals.  Downstream  of  the  bend,  however,  the 
rate  of  recovery  of  the  flow  is  too  slow,  which  points  to  the  same  weakness  in 
the  tubulence  model  as  found  in  the  recovery  region  of  the  flow  over  a  backward- 
facing  step. 
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NOMENCLATURE 


D 

De  (s  Re(0/Rc)1/2) 

k 

r 

Rc 

Re  (=  DUb/v) 

u*  (=  /  u2) 

u2 

U 


pipe  diameter 
Dean  number 

turbulent  kinetic  energy 
radial  (cross-stream)  coordinate 
curved  pipe  radius  of  curvature 
Reynolds  number 

rms  of  the  longitudinal  velocity  fluctuations 
longitudinal  component  of  normal  stress 
longitudinal  component  of  mean  velocity 
bulk  average  velocity  through  the  pipe 


v'  (s  /  vL 
,2 


vw 

V 

W 

X 


rms  of  the  ci rcumferenti al  velocity  fluctuations 

circumferential  component  of  normal  stress 

turbulent  shear  stress  in  r ,<t>  plane 

circumferential  component  of  mean  velocity 

radial  component  of  mean  velocity 

longitudinal  (streamwise)  coordinate  in  a  straight 
pipe  section 


Greek  letters 

e 

9 

v 

<t> 


dissipation  of  turbulent  kinetic  energy 

longitudinal  (streamwise)  coordinate  in  a  curved  pipe 
section 

fluid  kinematic  viscosity 
circumferential  (cross-stream)  coordinate 
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1.  INTRODUCTION 

1.1  The  Problem  of  Interest  and  Earlier  Work 

Over  the  last  ten  years,  approximately,  a  substantial  amount  of 
experimental  information  has  been  obtained  for  developing  turbulent  flows 
in  curved  ducts  of  rectangular  cross-section  [1-7].  Detailed  measure¬ 
ments  in  these  important  flow  configurations,  in  particular  of  the  secon¬ 
dary  motions  and  complex  turbulent  characteristics  which  they  induce, 
have  been  made  possible  by  the  non-intrusive  laser-Doppler  velocimetry 
technique.  In  contrast,  turbulent  flow-field  measurements  in  curved 
pipes  appear  to  be  relatively  sparse  [8],  This  is  surprising  given  the 
practical  importance  of  the  curved  pipe  configuration.  The  lack  of 
detailed  data  is  partly  explained  by  the  greater  difficulties  associated 
with  optical  alignment  through  transparent  curved  surfaces. 

The  dominant  influence  of  the  flat  walls  on  the  secondary  motion 
induced  in  rectangular  ducts  precludes  the  use  of  data  obtained  in  these 
configurations  to  interpret  accurately  the  characteristics  of  flow  deve¬ 
loping  in  curved  pipes.  For  example,  Cuming  [9]  shows  analytically 
that,  for  equal  radius  of  curvature  to  diameter  ratio  and  dimensionless 
streamwise  pressure  drop,  the  ratio  of  the  relative  intensity  of  the 
secondary  flow  in  a  duct  of  square  cross-section  to  that  in  a  pipe  of 
circular  cross-section  is  2.47  at  the  center  line  location.  This  result 
is  strictly  valid  only  for  fully  developed  laminar  flow.  Comparisons 
between  present  results  and  corresponding  measurements  by  Chang,  et  al . 
[5]  in  a  curved  duct  of  square  cross-section  with  similar  developing  tur¬ 
bulent  flow  characteristics  show  that  the  secondary  motion  ratio  can  be 
considerably  larger  than  2.47,  depending  on  the  streamwise  location. 

This  point  has  been  quantified  and  discussed  in  [10]. 


2. 


Among  the  early  measurements  of  secondary  turbulent  flow  patterns 
developing  in  and  downstream  of  curved  pipes  are  the  total  pressure  and 
yaw  results  obtained  by  Rowe  [11]  in  a  1B0°  bend  and  the  attached 
downstream  tangent.  In  Rowe's  experiment  Rc/D  =  12  and  Re  =  2.36  x  10^, 
where  Rc  is  the  bend  mean  radius  of  curvature,  D  its  diameter  and  Re  is 
the  Reynolds  number  based  on  the  bulk  fluid  velocity.  An  upstream 
tangent  of  length  X/D  =  69  was  used  to  provide  a  fully  developed  axi sym¬ 
metric  turbulent  velocity  profile  at  the  bend  inlet  plane.  Among  Rowe's 
findings  were  that:  secondary  flow  is  most  intense  at  about  0  =  30°,  as 
of  which  point  the  total  pressure  gradient  induces  a  streamwise  component 
of  vorticity  opposite  in  sense  of  rotation  to  the  streamwise  vorticity 
produced  at  the  start  of  the  bend;  the  curved  flow  is  essentially  fully 
developed  past  0  =90°;  there  is  evidence  of  local  reversal  in  the  secon¬ 
dary  flow  direction  along  the  bend  symmetry  plane  between  0  =  90°  and 
X/D  =  5. 

Because  of  the  nature  of  the  variables  measured  in  Rowe's  study,  and 
because  of  the  influence  of  the  mechanical  probes  used  in  the  flow,  an 
accurate  determination  of  the  fully  developed  flow  streamwise  location 
was  not  possible.  (A  combination  of  weak  secondary  motions,  undetected 
by  the  yawmeter,  with  unchanging  relatively  large  total  pressure  measure¬ 
ments,  would  certainly  have  given  the  impression  of  fully  developed 
flow.)  Due  to  experimental  uncertainty,  Rowe  was  not  able  to  define  pre¬ 
cisely  the  magnitude  and  extent  of  secondary  flow  reversal  in  the  vici¬ 
nity  of  the  bend  symmetry  plane.  Although  valuable  for  helping  to 
understand  the  behavior  of  the  mean  flow,  particularly  in  relation  to 
inviscid  effects,  the  pressure  and  yaw  contours  obtained  by  Rowe  shed  no 
light  on  the  turbulent  characteristics  of  the  flow. 


3. 


Enayet,  et  al  .  [12]  have  investigated  developing  turbulent  flow  in  a 
9(1°  curved  pipe  with  Rc/D  =  2.8  and  Re  =  43,000.  Relatively  thin  boun¬ 
dary  layers  were  induced  at  the  bend  inlet  plane  by  means  of  a  smooth 
contraction.  The  mean  longitudinal  (streamwise)  velocity  component  and 
the  corresponding  turbulence  intensity  were  measured  using  a  laser- 
Ooppler  velocimeter.  This  information  was  supplemented  by  wall-presure 
measurements.  The  relatively  flat  inlet  velocity  profile  and  small  value 
of  Rc/D  in  [12]  enhanced  streamwise  acceleration  of  the  flow  at  the  inner 
radius  wall  relative  to  the  outer  wall.  As  a  result,  throughout  the  bend 
and  up  to  X/D  =  6  in  the  downstream  tangent ,  the  maximum  in  the  stream- 
wise  velocity  component  remained  between  the  bend  center  line  and  the 
inner  radius  wall.  This  finding  contrasts  strongly  with  Rowe's  obser¬ 
vation  that  the  maximum  was  always  located  between  the  bend  center  line 
and  the  outer  radius  wall.  Because  the  Dean  numbers  of  these  two  flows 
are  very  similar,  being  De  =  Re  (D/Rc)1^2  =  6.8  x  104  in  Rowe's  bend  and 
2.6  x  104  in  the  bend  of  Enayet,  et  al.,  the  difference  in  the  maximum 
velocity  location  must  be  due  to  the  different  inlet  flow  boundary  layer 
thicknesses.  A  similar  finding  has  been  reported  in  [3]  for  the  flow 
through  a  curved  duct  of  square  cross  section. 

The  problem  of  numerical  simulation  of  developing  flow  in  curved  pipes  has 
been  the  subject  of  considerable  research  effort  for  the  past  twenty  years.  The 
extensive  literature  on  laminar  flow  has  been  recently  reviewed  in  [8,13],  The 
first  serious  attempt  to  compute  turbulent  flow  was  reported  by  Patankar  et  al . 
[14]  who  applied  their  three-dimensional  parabolic  (non-iterative)  procedure  to 
simulate  Rowe's  experiment.  The  parabolic  approximation  enforces,  at  any  sta¬ 
tion,  the  same  streamwise  pressure  over  the  whole  cross  section.  In  view  of  the 
relatively  large  radius  ratio  for  the  case  in  question,  such  a  numerical 
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approximation  was  satisfactory  over  much  of  the  bend.  At  the  inlet  to  and  exit 
from  the  bend,  however,  the  sign  of  the  pressure  gradient  is  in  fact  different 
on  the  inside  and  outside  of  the  bend  and  the  adjustments  to  the  flow  that  this 
causes  cannot  be  mimicked  by  a  parabolic  solver. 

Although  for  laminar  flow  several  workers  have  discretized  the  full  ellip¬ 
tic  form  of  the  equations  of  motion  f.15,16],  provided  there  is  no  flow  reversal 
in  the  streamwise  direction,  a  semi-elliptic  treatment  [17]  is  more  attractive 
as  only  the  pressure  field  then  requires  storing  over  the  solution  domain.  The 
approach  was  applied  by  its  originators  [17]  to  flow  in  a  rectangul ar-sectioned 
bend.  Recently  two  of  the  authors  [18]  have  applied  a  procedure  of  the  same 
type  to  predict  turbulent  flow  in  coils  and  in  the  90°  bend  of  [12].  For  the 
former  case  only  heat-transfer  data  were  available  for  comparison;  these  indi¬ 
cated  that  the  predictions  correctly  reproduced  the  monotonic  increase  in  local 
Nusselt  number  from  the  inside  to  the  outside  of  the  bend,  though  the  inner-to- 
outer  variation  was  only  about  60%  of  that  measured.  Agreement  with  the  stream- 
wise  velocity  profiles  of  [12]  was  broadly  satisfactory,  though  the  boundary 
layer  development  on  the  inside  of  the  bend  did  not  proceed  as  far  towards 
separation  as  the  experiment  showed,  nor  was  the  rate  of  recovery  downstream  as 
fast.  There  was,  moreover,  some  uncertainty  on  the  degree  of  difficulty  pro¬ 
vided  by  this  test  case  as  the  boundary  layers  were  rather  thin  and  the  pressure 
disturbance  so  large  that  inviscid  effects  would  have  dominated  over  much  of  the 
flow. 


1.2  The  Present  Contribution 

The  above  brief  review  has  Identified  substantial  deficiencies  in  our 
knowledge  of  the  physics  of  turbulent  flow  through  curved  pipes  and  of  the 
relaxation  processes  undergone  by  fluid  passing  from  a  curved  to  a  straight  pipe 
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section.  A  full  understanding  will  require  many  experiments  with  different 
ratios  of  bend  to  pipe  radii  and  different  flow  conditions  at  entry.  In  com¬ 
parison  with  [18],  we  have  selected  a  somewhat  larger  bend  radius,  180°  rather 
than  90°  of  turning,  and  a  fully-developed  flow  at  the  bend  entry  plane.  All 
these  differences  are  designed  to  make  the  resultant  flow  pattern  in  the  bend 
more  strongly  dependent  on  the  turbulent  stress  field.  Downstream  of  the  bend, 
the  turbulent  stresses  provide  the  only  substantial  agency  for  restoring  the 
flow  towards  that  associated  with  a  straight  pipe.  The  initial  stages  of  this 
restoration  have  been  documented  at  one  diameter  intervals. 

Finally,  in  order  to  gain  an  impression  of  how  successful  or  otherwise  the 
widely  used  k-e  eddy  viscosity  model  is  in  predicting  this  flow's  behaviour,  the 
solving  scheme  developed  in  [18]  has  provided  the  basis  for  a  numerical  simula¬ 
tion. 


2.  THE  EXPERIMENT 

2.1  Flow  System  and  Instrumentation 

The  experimental  system  was  composed  of:  a  water  rig  (described  in  detail 
in  [5]),  of  which  the  most  important  component  was  the  flow  test  section;  a  DISA 
55X  modular  series  laser-Doppler  velocimeter  and  its  associated  electronic 
instrumentation;  and  a  POP  11/34  Digital  Equipment  Corporation  minicomputer. 

The  basic  components  of  the  flow  test  section  are  shown  schematically  in 
Fig.  1.  They  comprised  two  straight  pipes  and  a  180°  curved  pipe,  constructed 
from  transparent  plexiglass.  The  pipe  cross-section  was  circular  throughout 
with  a  4,45  cm  +  0.02  cm  inner  diameter  (D).  The  ratio  of  bend  mean  radius  of 
curvature  to  pipe  diameter  was  Rc/D  =  3.375.  Both  tangents  were  of  length 
X  =  54.7  D,  being  respectively  attached  to  the  0°  (inlet)  and  180°  (outlet)  pla¬ 
nes  of  the  bend  by  means  of  flanges.  Special  care  was  taken  to  avoid  possible 
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mismatches  between  the  component  cross-sections  that  otherwise  might  disturb  the 
flow.  A  stainless  steel  (20  mesh)  screen  and  a  3.5  cm  long  flow  straightening 
honeycomb  section  were  placed  upstream  of  the  entrance  tangent.  They  were  held 
in  place  by  a  thin  circular  plexiglass  plate.  The  plate  was  3.175  mm  thick  and 
had  85  holes  of  3.175  mm  diameter  arranged  in  a  rectangular  array  spaced 
4.495  tun  on  the  centers  in  each  direction.  The  purpose  of  the  screen-honeycomb- 
plate  combination  was  to  make  uniform  the  cross-stream  plane  distribution  of  the 
flow  and  accelerate  its  streamwise  development. 

The  curved  pipe  section  was  constructed  by  fitting  together  two  symmetrical 
half  sections  of  plexiglass,  each  respectively  machined  on  one  of  the  two  flat 
faces  to  contain  the  shape  of  a  semicircular  open  channel.  This  method  of 
construction  ensured  that  when  matched  at  the  common  symmetry  plane  the  cross- 
section  of  the  resulting  curved  pipe  was  accurately  circular.  When  assembled, 
the  curved  pipe  section  had  flat  outer  surfaces  of  thickness  d  =  1.48  ±  0.02  cm 
(see  Fig.  1).  The  straight  pipe  sections  had  a  ci rcuwferenti ally  uniform  wall 
thickness  of  0.32  ±  0.02  cm.  The  entire  test  section  was  supported  horizontally 
by  a  metal  frame  with  the  bend  symmetry  plane  parallel  to  the  floor. 

The  test  section  was  part  of  a  closed  loop  through  which  water  was  made  to 
flow  by  gravity  from  a  constant  head  tank.  From  this  tank  the  flow  passed 
through  the  test  section,  and  then  into  a  large  sump  tank  from  where  it  was 
pumped  back  to  the  constant  head  tank.  Flow  to  the  head  tank  was  controlled  by 
a  gate  valve  and  measured  using  a  venturi  meter  connected  to  a  differential  mer¬ 
cury  manometer.  The  possibility  of  propagating  disturbances  through  the  test 
section  was  eliminated  by  avoiding  the  use  of  valves,  sharp  bends  and  metering 
devices  in  the  test  section  flow  loop.  Flow  to  and  from  the  test  section 
tangents  was  channeled  through  2  in  (5.08  cm)  i.d.  gently  bent  tygon  tube 
pieces.  Baffles  and  screens  located  in  the  constant  head  tank  served  to  dampen 
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the  swirling  motion  of  the  flow  leaving  the  tank.  Residual  swirl  and  weak 
secondary  motions  induced  by  the  mild  curvature  of  the  tygon  tube  upstream  of 
the  test  section  were  eliminated  by  the  flow  straightening  section. 

Measurements  of  the  mean  flow  and  turbulence  characteristics  were  made 
using  a  illSA  55X  modular  series  laser-Doppler  velocimeter  in  backscatter  mode. 
The  system  has  already  been  described  in  full  detail  in  reference  [5],  A  2-watt 
Lexel  Argon-Ion  water-cooled  laser  was  used  as  the  light  source.  The  laser  and 
velocimeter  optics  were  mounted  to  the  top  of  a  thick  aluminum  table  which  was 
itself  firmly  bolted  to  an  x,  y,  z  traversing  mechanism.  The  traversing  mecha¬ 
nism  could  displace  the  table  top  ±  7.5  cm  in  5  urn  increments  along  any  of  the 
coordinate  axes  by  means  of  three  linearly  encoded  stepping  motors  monitored  by 
the  POP  11/34  minicomputer.  The  minicomputer  functioned  as  the  central  data 
acquisition  and  reduction  controller  of  data  validated  and  measured  by  a  DISA 
55L96  Doppler  signal  procesor  or  "counter."  Directional  ambiguity  in  the  cir¬ 
cumferential  (or  tangential)  velocity  component  was  resolved  by  means  of  a  Bragg 
cell  combined  with  electronic  downmixing. 

2.2  Experimental  Methodology 

Prior  to  an  experimental  run,  water  was  allowed  to  flow  through  the  rig 
until  it  was  purged  of  air  bubbles  and  had  attained  a  steady  thermal  state 
corresponding  to  about  20°C.  Experiments  were  performed  for  two  values  of  the 
bulk  average  velocity  corresponding  to  Ub  =  1.29  ±  0.03  m/s  and 
Ufa  =  2.47  t  0.04  m/s,  respectively.  During  the  course  of  an  experimental  run 
the  fluid  temperature  rose  from  18°C  to  22°C,  approximately.  Values  for  the 
Reynolds  numbers  of  these  two  flows,  based  on  physical  properties  of  water  at 
20°C,  are  Re  =  57,400  and  110,000,  respectively.  The  associated  Dean  numbers 
are  De  =  31  ,300  and  59,900. 
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The  velocimeter  optical  probe  volume  was  formed  by  the  intersection  of  two 
514.5  nm  (green)  light  beams  with  a  half-angle  in  air  of  3.28°.  The  measurement 
volume  charateri sties  were:  a  diameter  of  0.1  mm,  a  length  of  1.6  mm,  and  a 
fringe  spacing  of  4.50  urn  with  about  18  fringes  contained  in  the  optical  probe. 
In  reality  spatial  filtering  and  threshold  settings  on  the  counter  reduced  the 
dimensions  of  the  optical  probe.  The  probe  volume  was  positioned  by  fine 
control  of  the  motorized  traversing  table.  The  details  of  the  optical  alignment 
procedure  are  provided  in  [10].  Cornstarch  particles  ranging  in  size  from  1  to 
10  urn  were  added  to  the  flow  so  that  validated  particle  data  rates  of  1  kHz  were 
achieved . 

With  a  reference  position  established,  usually  on  the  outer  surface  of  the 
test  section,  the  computer  software  was  activated  which  controlled  the  signal 
acquisition  and  data  processing.  An  initial  check  on  flow  symmetry  was  per¬ 
formed  at  X/D  =  -2  in  the  upstream  tangent.  Measurements  of  the  streamwise 
velocity  component  and  its  normal  stress  at  this  station  (half  profiles  are 
shown  in  Figs.  2  and  3)  are  in  good  agreement  with  the  data  of  laufer  [19], 
Measurements  of  the  circumferential  velocity  component  (not  plotted)  showed  that 
it  was  everywhere  less  than  ~  1.5»  of  the  bulk  average  velocity  at  this  station. 
Calculations  of  U^,  obtained  by  integrating  the  measured  velocity  distributions 
at  X/D  =  -2,  yielded  values  in  good  agreement  (to  within  ~  3%)  with  the  values 
obtained  from  the  venturi  meter  pressure  drop  readings. 

Following  the  symmetry  and  mass  flow  confirmations,  all  subsequent  measure¬ 
ments  were  restricted  to  vertical  (radial)  scans  in  the  surface  4>  =  jt/2  in  the 
symmetrical  upper  half  of  the  test  section.  Scans  were  made  at  X/D  =  -2,  -1,  1, 
2,  3,  4  and  5  in  the  straight  pipes  and  6  =  3°,  45°,  90°,  135°,  and  177°  in  the 
bend.  At  each  of  these  stations  15  radial  positions  were  probed,  starting  at 
the  centerline  and  moving  in  increments  of  1.5  mm  toward  the  wall,  with  a  16th 
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position  fixed  at  0.7  mm  from  the  inside  pipe  wall.  At  each  measurement  point 
the  mean  velocity  and  normal  stress  were  statistically  determined  from  popula¬ 
tions  of  1  to  3  samples  consisting  of  1,000  individual  measurements  each.  More 
than  90%  of  all  the  measurements  consisted  of  at  least  2  samples  per  point. 

Each  individual  measurement  was  required  to  satisfy  the  counter  5/8  validation 
comparison  to  within  a  preset  tolerance  of  6%.  At  every  validation  of  a  Doppler 
burst  a  "data  ready"  signal  was  issued  by  the  counter  to  a  logic  conversion  cir¬ 
cuit.  This  circuit  then  sent  a  triggering  pulse  to  the  computer  parallel  line 
interface  module  which  was  checked  for  data  availability  by  a  software  loop 
approximately  every  20  us. 

Although  the  velocimeter  is  capable  of  two-component  measurements,  cur¬ 
vature  of  the  test  section  surfaces  precluded  dual-channel  operation.  Values 
for  the  streamwise  component  of  the  mean  and  fluctuating  velocities,  U  and  u' , 
were  derived  from  measurements  obtained  with  the  probe  volume  fringes  aligned 
normal  to  the  longitudinal  coordinate  direction.  Values  for  the  circumferential 
components,  V  and  v' ,  were  obtained  through  a  rotation  of  the  optics  by  90°, 
thus  placing  the  fringes  normal  to  the  ci rcumferenti al  component  direction. 

Because  of  the  large  positive  magnitudes  of  the  streamwise  velocity  com¬ 
ponent,  frequency  shifting  was  unnecessary  for  measurements  of  this  component. 

By  contrast,  the  circumferential  velocity  component  exhibited  low  magnitudes, 
with  signs  both  positive  and  negative.  To  resolve  this  directional  ambiguity,  a 
net  frequency  shift  of  500  kHz  was  imposed  using  a  DISA  55N10  Bragg  cell  com¬ 
bined  with  electronic  downmixing. 

2.3  Measurement  error 

Error  sources  affecting  the  accuracy  (systematic  error)  and  precision 
(random  error)  of  1 aser-Dopp! er  measurements  have  been  discussed  in  [20-22].  In 
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this  study  the  most  serious  potential  systematic  errors  were  attributed  to  probe 
volume  positioning,  velocity  gradient  broadening  and  velocity  bias,  respec¬ 
tively.  The  main  sources  of  random  error  affecting  the  precision  of  the 
measurements  were  due  to:  a)  statistical  sampling  uncertainty  (due  to  the  finite 
size  of  sample  populations);  and  b)  uncertainty  in  the  determination  of  the  nor¬ 
malizing  bulk  average  velocity,  U^.  Estimates  of  the  first  uncertainty  were 
derived  from  the  measurements  themselves.  For  both  mean  velocity  components, 
this  uncertainty  was  always  less  than  ±  2-3%  (rms  error),  while  for  the  tur¬ 
bulence  intensities  it  was  always  less  than  t  2-5%.  The  error  in  Ub  was  t  2% 
and  arose  principally  from  uncertainities  in  the  manometer  pressure  readings. 

The  initial  probe  volume  location  uncertainty  was  estimated  to  be  half  the 
effective  probe  volume  length,  or  ±  0.5  irni  approximately.  This  error  is  altered 
as  the  beams  traverse  the  pipe  because  of  changes  in  index  of  refraction  (at  the 
air/plexiglass  and  plexiglass/water  interfaces)  and  wall  curvature  effects.  To 
calculate  the  actual  probe  volume  location,  and  to  estimate  the  uncertainty  in 
this  location,  it  is  necessary  to  know  the  relevant  trigonometric  relations. 
These  are  derived  in  [10].  The  relations  depend  on  the  velocity  component  of 
interest  and  on  the  shape  of  the  surfaces  at  the  measurement  station  in 
question.  In  regions  of  surface  curvature,  both  the  probe  volume  location  and 
the  half  angle  between  the  beams  in  water  depend  non-linearly  on  the  probe 
volume  virtual  radial  location  (the  location  in  air).  Calculations  of  the  error 
incurred  in  y,  the  half  angle  of  the  beams  in  water,  due  to  a  ±  0.5  mm  posi¬ 
tioning  error  in  the  initial  probe  volume  location  showed  that  this  source  of 
inaccuracy  was  negligible.  Similarly,  calculations  showed  that  the  initial 
positioning  uncertainty  was  only  significant  very  near  the  pipe  walls,  in 
regions  with  steep  field  gradients.  Thus,  estimates  of  the  maximum  (absolute) 
systematic  errors  in  the  measurements  obtained  nearest  the  walls  gave:  ~  4%  for 
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U/Uh;  ~  10%  for  u7Ub;  ~  15%  for  V/Ub ,  and  ~  10%  for  v’/Ub.  In  the  bulk  of  the 
flow  all  these  errors  were  significantly  less  than  ~  2%. 

Velocity  gradient  broadening  has  been  analyzed  by  Mel  ling  [23]  who  proposed 
a  simple  method  for  estimating  its  magnitude.  For  the  conditions  of  this  study, 
gradient  broadening  in  the  measurements  was  estimated  to  be  less  than  1% 
everywhere  in  the  flow. 

Various  weighting  methods  have  been  proposed  to  correct  for  the  velocity 
bias  effect  [21,22],  but  none  of  these  is  entirely  satisfactory.  They  all 
involve  assumptions  regarding  the  statistical  distribution  of  particles  in  the 
flow.  In  addition,  the  corrections  can  be  influenced  by  the  problem  of 
"incomplete-signal  bias."  The  use  of  the  weighting  scheme  proposed  by 
McLaughlin  and  Tiederman  [24],  in  which  the  particle  residence  time  is  taken  to 
be  inversely  proportional  to  the  modulus  of  the  measured  velocity,  was  investi¬ 
gated  in  this  study.  For  the  longitudinal  velocity  component,  differences  bet¬ 
ween  weighted  and  unweighted  data  sets  were  relatively  small  (less  than  1%). 

For  the  ci rcumferential  component  the  differences  were  much  larger,  with  the 
weighted  data  giving  a  less  realistic  (overdamped)  picture  of  the  flow  field. 

The  reason  for  this  is  attributed  to  an  overcorrection  in  the  weighting  scheme 
which,  for  the  circumferential  component,  amounts  to  a  di sproportionate 
weighting  of  slow  moving  particles.  This  occurs  as  a  consequence  of  taking  the 
particle  residence  time  as  t  „  ~  V'*  instead  of  t  ~  (U2  +  V2  +  W2)'1,/2.  The 
latter,  more  accurate,  weighting  could  not  be  performed  with  the  present  single 
component  measurement  system. 

For  turbulence  intensities  less  than  10%  Drain  [21]  shows  that  velocity 
bias  corrections  of  the  mean  and  r.m.s.  velocities  should  not  amount  to  more 
than  ~  2%.  In  this  study  the  measured  turbulence  intensities  never  exceed 
~  12%.  As  a  result,  plots  of  the  unweighted  data  only  are  presented  here. 
Listings  of  both  weighted  and  unweighted  measurements  are  provided  in  [10]. 
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3.  THE  COMPUTATIONAL  METHOD 

The  numerical  procedure  developed  by  the  authors  to  predict  laminar  or  tur¬ 
bulent  flow  in  pipe  bends  has  been  described  in  [13,18].  Here,  therefore,  we 
simply  summarize  its  main  features. 

The  method  is  based  on  a  semi-el  1 iptic  discretization  of  the  Reynolds 
equations  expressed  in  toroidal  coordinates  in  which  the  three  components  of 
velocity  and  the  static  pressure  form  the  set  of  dependent  variables.  Following 
the  methodology  of  the  TEACH  family  of  computer  programs,  [25],  a  finite-volume 
approach  to  forming  the  difference  equations  is  adopted  with  a  staggered  mesh 
for  the  different  dependent  variables.  The  quadratic  upstream  weighted  approxi¬ 
mation  is  adopted  for  convective  transport,  [26],  while  a  modified  version  of 
Patankar's  SIMPLER  algorithm  [27]  is  used  to  secure  a  pressure  field  in 
compliance  with  the  continuity  equation. 

Although  in  a  shear  flow  as  complex  as  this  one  might  well  expect  to  have 
to  adopt  some  non-equilibrium,  non-isotropic  model  of  turbulence  to  secure  close 
agreement  with  experiment,  the  (relative)  computational  simplicity  of  using 
isotropic  effective  diffusion  coefficients  persuaded  us  that  in  this  study  the 
turbulent  stresses  should  be  obtained  by  the  standard  k~e  model  [18,28], 

However,  in  place  of  the  wall  functions  that  usually  provide  a  link  between 
near-wall  velocities  and  shear  stresses  and  supply  effective  boundary  conditions 
for  the  turbulence  energy  (k)  and  its  viscous  dissipation  (t ) ,  a  fine  mesh  has 
been  incorporated  to  cover  the  viscous  sublayer  and  "buffer"  regions.  In  this 
sub-region  we  have  for  simplicity  adopted  the  mixing-length  hypothesis  using  van 
Oriest's  [29]  damping  function.  Of  course,  to  resolve  accurately  the  very  rapid 
variation  of  effective  transport  coefficients  across  this  near-wall  sublayer  a 
fine  mesh  is  required.  In  order  that  this  requirement  should  not  impose  a 
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substantial  core  penalty,  the  variation  of  pressure  across  the  sublayer  was 
obtained  by  radial  equilibrium  (instead  of  by  continuity);  the  technique,  known 
as  PSL,  is  described  in  [30]. 

Computations  were  made  over  the  semi-circular  cross  section  bounded  by  the 
plane  of  symmetry  passing  through  the  center  of  the  pipe  and  the  bend  center. 
Starting  inlet  profiles  were  obtained  by  computing  the  development  of  the  flow 
in  a  straight  pipe  over  a  distance  of  53  diameters.  The  downstream  profiles  of 
the  dependent  variables  were  then  used  as  the  inlet  profiles  for  the  computation 
of  the  U-bend  flow  at  a  position  2  diameters  upstream  of  the  bend.  The  mesh 
comprised  28  radial  and  24  uni fcrmly-spaced  ci rcumferenti al  nodes  mapping  the 
semi-circular  section  and  160  streamwise  planes.  This  mesh  density  has  been 
found  from  our  previous  work  to  give  essentially  grid  independent  resultsl  when 
used  with  the  third-order  accurate  quadratic  upstream  approximation  for  convec¬ 
tion.  The  computations  extended  a  total  of  8  diameters  into  the  tangent 
downstream  of  the  bend  where,  at  the  exit  plane,  zero-streamwise  gradient 
constraints  were  applied  to  all  variables.  Successive  streamwise  passes  over 
the  domain  were  made  until  momentum  and  mass  residuals  summed  over  all  planes 
were  below  1%  of  the  entering  momentum  and  mass  flow  rates.  Convergence 
required  about  6  hours  cp  time  on  a  Cyber  205  (with  an  unvectorized  code) 
though  in  practice  only  very  small  changes  resulted  over  the  second  half  of  the 
computati on . 


♦♦**★★***★****★****★★ 


lFor  example,  in  [18]  an  increase  in  the  number  of  radial  nodes  from  28 
to  40  altered  the  predicted  level  of  Nusselt  number  by  less  than  U. 
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4.  RESULTS  AND  DISCUSSION 

Measurements  of  the  longitudinal  (U/U^)  and  circumferential  ( V/U^)  velocity 
components  and  their  respective  turbulence  intensities  (u'/U^,  v'/Ub)  are  shown 
in  Figs.  2-3.  Only  the  data  at  the  lower  Reynolds  number  (Re  =  57,400)  are  here 
included.  Except  for  relatively  small  differences,  the  dimensionless  profiles 
show  that  flow  development  is  essentially  the  same  at  Re  =  110,000  [10].  This 
is  an  important  point  since  all  the  calculations  reported  here  were  performed 
for  Re  =  65,000  (with  Rc/D  =  3.375).  Tabulations  of  the  measurements  are 
available  in  [10]. 

The  U  and  u'  component  measurements  at  X/D  =  -2  show  good  agreement  with 
the  developed  flow  measurements  obtained  by  Laufer  [19]  in  a  straight  pipe.  At 
X/D  =  -2  the  V  component  is  less  than  (+)  0.015  Ub  (not  plotted).  However,  at 
X/D  =  -1  and  0=3°,  the  V  component  is  ~  (+)  0.04  U^.  With  reference  to  the 
coordinate  system  in  Fig.  1,  this  implies  that  the  inlet  flow  is  gradually  (but 
only  slightly)  accelerated  towards  the  inner  radius  of  the  pipe,  a  result 
expected  from  potential  flow  theory  and  previously  observed  by  others. 

In  the  first  half  of  the  bend  (0  =  3°  to  90°),  the  U-profiles  show  the  core 
of  the  streamwise  flow  losing  speed  while  the  flow  near  the  wall  accelerates. 
V-profiles  at  corresponding  stations  reveal  the  development  of  a  strong  secon¬ 
dary  flow.  This  secondary  flow  is  induced  by  the  transverse  pressure  gradient 
set  up  between  the  outer  (rQ)  and  inner  (r^)  wall  regions  of  the  bend.  In  the 
pipe  center,  it  works  to  overcome  and  reverse  the  sense  of  the  cross-stream 
motion  in  the  inlet  flow.  Near  the  wall  the  original  sense  of  cross-stream 
motion  is  preserved  and  intensified,  with  (V/Ub)max  ~  0.30  at  8  =  45°. 
Nevertheless,  the  levels  of  secondary  flow  found  in  the  present  study  are  lower 
than  those  obtained  in  a  square  sectioned  duct  for  equivalent  curvature  ratio 
and  entering  flow  conditions  [5], 
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Between  e  =  45°  and  135°  the  V-profiles  reveal  a  striking  feature  of  the 
flow.  In  the  fluid  core,  r/(D/2)  <  0.5,  the  cross-stream  flow  undergoes  a 
second  reversal  in  its  sense  of  motion  and  is  redirected  towards  the  pipe  inner 
radius.  As  a  result,  a  region  of  negative  cross-stream  flow  (directed  from 
r.  towards  r  )  is  trapped  between  the  core  and  the  wall;  for  example, 

V/Ilp  ~  -0.06  at  e  =  135°.  Between  8  =  135°  in  the  bend  and  X/D  =  5  in  the 
downstream  tangent,  the  region  of  negative  cross-stream  flow  gradually  disap¬ 
pears.  Notwithstanding,  an  imprint  of  its  presence  remains  in  all  the  sub¬ 
sequent  V-profiles. 

The  occurrence  of  a  second  cross-stream  flow  reversal  past  8  =  90°  in  the 
bend  supports  the  concept  proposed  by  Rowe  [11],  that  each  symmetrical  half  sec¬ 
tion  of  the  bend  develops  two  counter-rotating  vortical  structures.  The  more 
intense  of  the  two,  located  between  the  pipe  wall  and  the  core  of  the  flow,  pre¬ 
serves  the  sense  of  cross-stream  motion  induced  by  the  transverse  pressure  gra¬ 
dient  at  the  start  of  the  bend.  The  smaller,  weaker,  structure  is  mainly 
confined  to  the  core  and  is  attributed  to  the  formation  of  a  transverse  pressure 
gradient  opposite  in  sign  to  that  at  the  start  of  the  bend. 

The  present  mean  flow  measurements  place  on  a  firm  quantitative  basis  the 
interpretation  by  Rowe  of  his  total  pressure  and  yaw  measurements.  However,  in 
contrast  to  Rowe's  results,  the  V-profiles  provided  here  show  that  fully  deve¬ 
loped  mean  flow  conditions  are  not  established  in  the  present  bend. 

The  computed  distribution  of  longitudinal  velocity  shown  in  Fig.  2  actually 
mirrors  very  closely  the  measured  development.  Only  at  90°  is  there  a  discer¬ 
nible  difference  between  the  two,  the  experiments  indicating  a  larger  dip  in 
velocity  at  the  axis.  This  measure  of  agreement,  while  encouraging,  at  least 
partly  reflects  that  the  mean  velocity  profile  along  this  $  =  ir/2  surface  does 
not  offer  a  particularly  sensitive  test.  This  is  not  the  case  with  the  secon- 
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dary  flow  shown  in  Mg.  2.  The  computations  do  not  reproduce  as  strong  an 
increase  in  the  V  velocity  component,  directed  towards  the  inside  of  the  bend, 
as  the  data  at  X/D  =  -1  and  3  =  3°.  However  they  do  predict  very  well  the 
strong,  single-cell  secondary  flow  that  develops  by  45°.  They  also  capture 
broadly  the  way  the  "return"  flow  near  the  mid-plane  is  brought  to  a  halt  and 
reverses  by  135°.  Interestingly,  the  computations  suggest  that  the  sign  of  the 
secondary  velocity  on  the  symmetry  plane  reverses  yet  again  at  177°  and  at  one 
diameter  downstream.  While  the  measurements  at  177°  still  show  marginally  posi¬ 
tive  V  for  small  r,  the  sign  of  the  acceleration  is  evidently  correct.  The  dif¬ 
fering  development  of  the  secondary  flow  at  different  radial  positions  is 
brought  out  in  Fig.  4  which  shows  the  progress  of  V/U  on  the  n/2  surface  at 
2r/D  =  0.75  and  at  the  pipe  axis.  The  flow  reversals  exhibited  by  the  measure¬ 
ments- are  generally  well  captured  by  the  numerical  computation,  at  least  within 
the  bend  itself. 

A  more  complete  picture  of  the  computed  flow  development  in  the  U-bend  is 
provided  in  Fig.  5  which  shows  both  the  secondary  flow  vectors  in  the  cross- 
sectional  plane  and  the  streamwise  flow  contours.  The  modifications  that  the 
former  bring  to  the  latter  are  readily  discernible.  Particularly  interesting 
are  the  90°  a  '  135°  contours  where  the  maximum  streamwise  velocity  is  displaced 
off  axis  (as  in  the  experiment,  Fig.  2),  and  those  at  177°  where  the  contours 
U/Ub  =  0.8  and  1.0  show  clearly  the  effect  of  the  reversal  of  the  secondary  flow 
on  the  mid-plane.  In  Fig,  5,  R  and  L  are  respectively  the  bend  radius  of  cur¬ 
vature  and  length  of  the  pipe  "centerline"  measured  from  the  inlet  plane. 

Downstream  of  the  bend  the  computations  do  not  mimic  the  steady  cross- 
stream  movement  of  fluid  towards  the  side  of  the  pipe  corresponding  with  the 
inside  of  the  bend  (Fig.  3).  This  reflects  a  significantly  too  slow  recovery  of 
the  computed  streamwise  velocity  from  the  strong  adverse  pressure  gradient  that 
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the  fluid  on  the  inside  of  the  bend  encounters  at  the  bend  exit.  Computations 
for  the  90°  bend  of  Enayet  [18]  with  precisely  the  same  model  show  a  similar  too 
slow  recovery.  It  is  also  relevant  to  mention  that  several  groups  adopting  the 
same  (or  a  very  similar)  turbulence  model  to  compute  the  flow  downstream  of  the 
reattachment  point  behind  a  backward-facing  step  for  the  1981  Stanford 
Conference  also  found  a  too  sluggish  recovery  [31].  The  behaviour  here  pre¬ 
dicted  thus  seems  an  endemic  feature  of  the  model  that  needs  to  be  removed. 

In  flow  through  a  straight  pipe,  the  circumferential  normal  stress  is 
approximately  equal  to  the  arithmetic  mean  of  the  other  two  components. 
Accordingly,  on  Fig.  3,  showing  the  measured  distribution  of  v' ,  have  been  added 
the  profiles  of  /  2k/3  using  the  k-levels  generated  by  the  computations.  As 
expected,  at  the  inlet  there  is  close  agreement  between  computation  and  measure¬ 
ment.  There  is,  moreover,  reasonably  close  agreement  between  the  two  over  the 
whole  flow  domain.  We  note,  however,  that  measured  levels  of  v'  are  rather 
higher  than  those  computed  near  the  wall  between  45°  and  177°.  This  difference, 
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we  believe,  arises  mainly  from  the  direct  generation  of  vc  by  the  strongly 
sheared  secondary  flow  (i.e.,  large  3V/3r);  although  generation  of  k  by  this 
agency  also  appears  in  the  turbulence  energy  transport  equation,  it  is  a  less 
important  process  there.  At  90°  and  135°  the  computed  levels  of  ✓  2k/ 3  are 
lower  than  the  measured  levels  of  v'  over  the  whole  radius  by  amounts  ranging 
from  1  to  30%.  It  seems  likely  that  at  these  positions  the  computed  energy 
dissipation  rate,  e,  has  become  too  large,  preventing  the  rise  in  fluctuating 
velocity  levels  that  the  experimental  values  indicate.  The  experimental  values 
of  v'  decrease  markedly  beyond  the  bend  exit,  and  the  subsequent  downstream 
development  shows  remarkable  agreement  between  the  computations  and  the  measure¬ 
ments. 

Plots  of  the  ratio  v'/u',  a  measure  of  turbulent  flow  anisotropy,  are  given 
in  [10].  The  results  suggest  that,  between  X/D  =  -1  and  8  =  90°,  near  the  outer 


radius  wall,  the  term  -  vw  3V/a  r  contributes  strongly  to  the  generation  of  v2. 

Beyond  0  =  90°  this  contribution  decreases ,  allowing  dissipation  and  redistribu- 
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tion  mechanisms  to  reduce  the  level  of  v  near  the  wall,  as  observed. 

5.  CONCLUDING  REMARKS 

This  study  places  on  a  firm  quantitative  basis  some  of  Rowe’s  [11]  qualita¬ 
tive  interpretations  of  the  behaviour  of  turbulent  flow  through  a  180°  curved 
pipe.  Specifically,  the  existence  of  two  cross  stream  flow  reversals  has  been 
confirmed  both  experimentally  and  by  the  numerical  simulations.  The  results 
support  the  notion  of  an  additional  (symmetrical)  pair  of  counter-rotating  vor¬ 
tical  structures  embedded  in  the  core  of  the  flow  within  the  curved  pipe.  The 
levels  of  secondary  flow  are  on  the  whole  lower  than  in  the  corresponding  square 
duct  flow  [5].  This  underlines  the  dominant  effect  of  the  flat  walls  in  the 
curved  square  duct  in  driving  the  secondary  motion.  No  cross-stream  flow  rever¬ 
sals  were  measured  in  the  latter  case,  which  further  indicates  how  different  are 
the  structures  in  what  superficially  appear  to  be  two  very  similar  flows. 

The  rms  level  of  velocity  fluctuations  is  substantially  raised  as  the  flow 
passe.;  around  the  bend  due  to  the  additional  mean  strain  associated  with  the 
turning  of  the  primary  flow  and,  particularly,  the  secondary-flow  velocity  gra¬ 
dients  thus  created.  There  are,  in  addition,  changes  in  the  structure  of  the 
turbulent  field  in  the  core  -  notably  the  large  increase  of  streamwise  fluc¬ 
tuations  at  90°  and  the  rise  in  v'  observed  between  135  and  177°  -  whose  cause 
cannot  be  unravelled  with  any  certainty  from  measurements  on  the  $  =  n/2  surface 
al  one. 

The  numerical  simulations  have  reproduced  with  a  gratifying  degree  of  fide¬ 
lity  the  measured  evolution  of  the  flow.  It  strengthens  the  impression  formed 
in  [5]  that  often  the  k~e  eddy  viscosity  model  does  better  in  simulating  real  ly 
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complex  flows,  such  as  those  found  within  the  180°  bend,  than  it  does  in  less 
strongly  perturbed  flows.  For  example,  in  the  present  case  the  prediction  of 
the  flow  recovery  downstream  of  the  bend  shows  the  largest  disagreement  with  the 
measured  behavior  even  though  secondary  velocities  are  there  very  weak.  No 
general  conclusions  may  be  drawn,  however,  for  computations  of  the  square  bend 
with  the  same  model  of  turbulence  used  here  [32,33]  display  spectacular 
disagreement  with  the  data  of  [5]. 

Future  work  on  bend  flows  needs  to  be  aimed  at  providing  a  more  complete 
mapping  of  the  primary  and  secondary  mean  flow  and  stress  distributions  over  the 
entire  conduit  cross-secti on .  This  data  should  help  in  clarifying  the  role  of 
the  additional  pair  of  vortical  structures  in  curved  pipe  flow  and  the  relaxa¬ 
tion  process  occurring  in  the  downstream  tangent. 
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Fiq.  3  Measurements  (points)  and  calculations  (continuous  lines)  of  the 

longitudinal  (u ‘ )  and  circumferential  (v1)  turbulence  intensity  com¬ 
ponents  at  sequential  longitudinal  stations  in  a  180  curved  pipe  with 
straight  tangents.  Dashed  line  corresponds  to  Laufer  s  [19]  developed 
flow  data.  In  the  calculations  u‘  =  v'  ■  (2k/3) 1/2 . 


Fig.  4.  Measurements  (points)  and  calculations  (lines)  of  developing  secondary 
flow  at  two  radial  locations  contained  in  the  symmetry  plane,  $  = 
of  a  180°  curved  pipe  with  straight  tangents. 
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Proc.  4th  Int.  Conf.  on  Numerical  Methods  in 
Laminar  and  Turbulent  Flow,  pp. 1023-1045, 
Swansea,  July  9-12,  1985 


;  ASM  PREDICTIONS  OF  TURBULENT  MOMENTUM  AND  HEAT  TRANSFER  IN 
j COILS  AND  U-BENDS 
I  Author  Heie 

H.  Iacovides  and  B.E.  Launder 
Mechanical  Engineering  Department, 

! UMIST,  Manchester  M60  1QD 

I 

i 

I  first  Tiie  papet; describes  the  incorporation  of  an  algebraic 
stress  model  (ASM)  of  turbulence  into  an  accurate 
semi-elliptic  solution  procedure  for  the  Reynolds  equations 
and  the  corresponding  energy  eauation.  Particular  attention 
is  given  to  the  devices  necessary  to  procure  convergent 
behaviour.  Applications  of  the  approach  are  reported  to  the 
computation  of  fully-developed  flow  in  coils  and  of  flow 
developing  in  U-bends.  Generally  encouraging  agreement  with 
Ithe  experimental  data  is  obtained  though  the  computations  do 
not  mimic  the  partial  laminarization  that  appears  to  occur  in 
the  strong  acceleration  on  the  inside  of  the  bend  at  inlet. 

I.  INTRODUCTION 

j  - 

The  numerical  simulation  of  developing  turbulent  flow 
through  tightly  curved  pipes  poses  a  number  of  problems 
concerning  both  the  physics  -  the  model  for  turbulent 
transport  -  and  the  numerics.  The  flow  is  highly  three- 
dimensional,  for  the  centripetal  acceleration  of  the  sheared 
jstreamwise  flow  induces  a  strong  secondary  motion  with  peak 
(Velocities  30%  or  more  of  the  bulk  streamwise  velocity.  In 
fact,  the  secondary  flow  is  strongest  very  close  to  the  pipe 
walls  where  cross  currents  steadily  carry  fluid  from  the 
outside  to  the  inside  of  the  bend.  This  feature  means  that 
the  near-wall  velocity  profile  is  unlikely  to  conform  closely 
with  the  usual  semi-logarithmic  law.  Accordingly,  the  normal 
practice  of  making  the  velocity  field  comply  not  with  the 
no-slip  boundary  condition  at  the  wall  but  with  the  semi- 
logarithmic  profile  at  some  appropriate  distance  from  it 
needs  to  be  abandoned. 

This  need  accentuates  the  problem  of  sheer  numerical 
resolution.  For  example,  to  compute  the  turbulent  flow 
around  a  U-bend  vhtltout  substantial  numerical  error  requires, 
with  3rd-order  quadratic  upwind  differencing  and  a  fine  near- 
wall  mesh  to  resolve  the  viscous  sublayer,  the  velocities- to 
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be  computed  at  some  50  thousand  points  in  the  domain.  With 
such  a  large  number  of  mesh  points  it  is  not  routinely 
possible  to  store  in  core  all  the  dependent  variables 
'  (pressure,  three  velocity  components  and  usually  two 
{turbulence-characterizing  parameters).  Fortunately,  if  t 
; streamwice  flow  recirculation  is  absent,  a  semi-elliptic 
treatment  may  be  adopted  in  which  only  the  pressure  is  stored 
j  over  the  whole  domain  while  the  values  of  all  other  variables 
;  are  held  just  on  two  adjacent  cross-sectional  planes.  This 
i  approach,  developed  by  Pratap  and  Spalding  [1],  has  been 
applied  to  flow  around  square  bends  by  its  originators  [1], 
j  Chang  et  al  [2] ,  and  Birch  [3]  and  to  laminar  and  turbulent 
,  flow  through  circular-sectioned  pipes  by  the  authors  [4,5,6]. 

Author-  The  computational  results  for  turbulent  flow  presented  1 
in  ref [5]  were  based  on  a  two-region  turbulence  model.  Over  i 
most  of  the  flow  the  widely  used  k~ e  eddy  viscosity  model 
|  (EVM)  was  adopted,  the  value  of  the  turbulent  viscosity  at 
any  point' being  obtained  from  local  values  of  the  turbulence  ; 
i  energy  (k)  and  its  dissipation  rate  (e).  Very  near  the  wall, 

1  however,  where  viscous  effects  become  of  direct  importance. 
Van  Driest's  version  of  the  mixing-length  hypothesis  was 
used.  The  resultant  computations  displayed  fairly 
satisfactory  agreement  with  experiment  although  systematic. 

;  differences  were  nevertheless  evident.  Predictions  of  local 
■  heat  transfer  coefficient,  while  giving  a  mean  level  broadly 
in  line  with  experiment,  exhibited  a  too  small  variation 
I  between  the  inside  and  outside  of  the  bend.  Qualitatively  , 

!  this  shortcoming  was  foreseen.  Eddy  viscosity  schemes  are 
I  known  not  to  capture  correctly  the  strong  damping  or 
,  augmentation  of  turbulence  associated  with  streamline 
curvature  [7] . 

Now,  in  simpler,  two-dimensional  strains  it  is  well 
I  established  that  turbulence  models  based  on  the  Reynolds 
stress  transport  equations  do  succeed  fairly  well  in  mimick- 
i  ing  the  observed  strong  effects  of  curved  streamlines  in  the 
|  turbulent  stresses.  Our  objective  in  the  work  reported  ! 

|  hereunder  has  thus  been  to  incorporate  one  of  these  higher  I 
level  turbulence  models  into  the  solving  procedure  developed 
in  ref  [5] .  We  use  what  is  known  as  an  "algebraic  stress 

model"  (ASM)  which  truncates  the  Reynolds  stress  transport 
,  equations  to  a  direct  algebraic  connection  between  the  local  , 

I  Reynolds  stresses  and  the  local  strain  field  -  though  that 
*  connection  is  far  more  intricate  than  with  eddy  viscosity  ! 

1  models. 

I 

j  This  enhancement  of  turbulence  model  has  necessitated 

j  major  changes  to  the  way  the  solution  procedure  is  organized 
in  order  to  secure  convergent  behaviour.  An  account  of  the 
novel  practices  introduced  is  given.  In  view  of  length  ' 
limits  on  the  pr.per,  aspects  of  the  procedure  which  remain 
'  .unchanged  from  our  earlier  work  [4,5,6]  are  _summari zed  only 
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briefly.  Comparisons  are  drawn  in  Section  4  with  the 
available  experimental  data  and  with  computations  generated 
,  with  the  k~e  EVM. 


j  2.  THE  MATHEMATICAL  AND  PHYSICAL  MODEL 
‘2.1  The  Describing  Mean  Flow  Equations 

I 

1  The  Reynolds  equations  describing  the  motion  of  an 

;  incompressible,  stationary,  turbulent  flow  in  toroidal 
•coordinates.  Figure  1,  may  be  expressed 

I  P  |  C(V)  +  Sc(f)j  -  D(40  -  p  J  R(Y)  +  sR(t')|  +  sd(T)  +  Sp(V 


(1) 

wher'^  frstands,'for°  any  of  the  components  of  mean  velocity  u 
I  (azimuthal,  $  ) ,  V  (radial,  r)  or  W  (streamwise,  e  ).  The 
quantities  C,  D  and  R  denote  the  convection,  molecular- 


diffusion  and  turbulent-diffusion  operators 
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where  V  refers  to  the  fluctuating  velocity  component  in  the 
I  same  direction  as  ¥  and  r  is  the  local  radius  from  the 
centre  of  the  bend  ( R  +  rccos  e )  ,  R  being  the  mean  bend 
radius.  (When  considering  flow  in  straight  tangent  sections 
;  to  a  bend  R  is  set  to  same  appropriately  large  value) .  i 

1  } 
i  The  remaining  source  terms  in  eq(l)  are  set  out  in  Table  j 
j  1. 

Equation  (1)  also  describes  the  mean  energy  equation  for  ! 
low-speed  flow  with  ?  now  standing  for  temperature,  T  (and 
the  fluctuating  temperature,  t) .  In  that  event,  Cff)  and 
1  R ( 'F)  retain  precisely  the  same  form  as  for  the  velocity 
i  components  vrfiile  in  D(T)  the  viscosity  v  is  replaced  by  v/o, 
j  a  denoting  the  fluid  Prandtl  number.  In  the  energy  equation  j 
1  the  residual  soirees  are  all  zero.  | 


j  2.2  The  Turbulence  Models 

(a)  The  high  Reynolds  number  region  :  From  the  pipe  axis  to 
a  radius  equal  to  96%  of  the  pipe  radius  one  of  two 
turbulence  models  are  adopted,  the  standard  k~  e  eddy 
.viscosity  model  (EW)  and^  an  algebraic  stress  model  (ASM) . 
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The  EVM  adopts  the  Newtonian  stress-strain  law 
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1  J 


3U£ 

3x- 


+  |SL} 

3x. 

1 


(2) 


I  (u.  and  U.  denoting  fluctuating  and  mean  velocities)  where 
the  turbulent  viscosity  is  expressed  in  terms  of  the 
;  turbulent  kinetic  energy  k  and  its  rate  of  dissipation  e: 

!  v  «  0.09  k2/e 

j  nrib  OF  ?v  !  *  r. '  <  r*  t;  R  r 

and  k  and  e  are  obtained  from  transport  equations  as 
described  below.  The  turbulent  heat  fluxes  u^t  are  likewise 
.  obtained  from 
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I  where  oT  ,  the  turbulent  Prandtl  number  for  energy  diffusion, 
j  }s  taken  as  ,0.9..,;  j 

The  ASM  adopted  in  this  study  is  essentially  the  usual 
isotropization-of-production  closure  with  Rodi's  [8] 
transport  hypothesis 
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1  In  equation  (4)  P. .  is  the  stress  generation  rate  by  mean 
.  strain  ** 
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and  4>ij2  are  wall-proximity  functions  accounting  for 
pressure  reflections  from  the  rigid  walls  of  the  pipe.  We 
adopt  the  forms  suggested  by  Shir  [9]  and  Gibson  and  Launder 
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n.  is  the  unit  vector  normal  to  the  wall,  x^  is  the  normal 
distance  from  the  wall  and  l  is  the  turbulent  length  scale 
0.38  k^2  /  e  (the  numerical  coefficient  being  chosen  so  that 
in  an  equilibrium  near-wall  layer  the  quantity  l/x^  is 
unity.  The  function  f  <l/x  )  was  simply  taken  as  i/xj!  by 
Gibson  and  Launder-.  - . '-*>  - -  - -> 


Reference  (10J  had  proposed  the  above  form  for  a 
turbulent  shear  flow  near  an  infinite  plane  surface  (the 
atmospheric  boundary  layer) .  The  problem  of  how  to  adapt  the 
approach  to  apply,  with  any  degree  of  rigour,  within  a 
toroidal  enclosure  is  a  difficult  one  that  we  have  not 
j  attempted  to  address  here.  We  have  for  simplicity  retained 
;  the  form  given  above  merely  replacing 

!  .x  r  .1  -fl 

!  x.  by  (a  -  r)  -  (a  +  r) 

•  riTLC  OF  PAI'LF  "  L  -1 

where  a  is  the  pipe  radius.  Close  to  the  wall  the  second 
term  in  square  brackets  makes  only  a  small  contribution 
i  compared  with  the  first  which  conforms  with  our  expectation 
1  that  turbulence  there  should  know  little  about  the  overall 
geometry  of  the  cross  section.  At  the  axis,  however,  the  two 
terms  cancel  thus  removing  any  wall  correction.  This  feature 

!  is . necessary  when  working  within  the  framework  of  the 

i  indicated  model'  for  otherwise  at  r  *  0  we  should  get 
different  values  for  u2  and  v2.  While  not  regarding  the 
;  above  simple  expedient  with  much  satisfaction,  it  is  our 
;  impression  that  its  weaknesses  are  not  having  a  great  effect 
on  the  calculated  flow  and  thermal  behaviour. 

To  calculate  the  turbulent  heat  fluxes  we  have  adopted 
the  generalized  gradient  diffusion  hypothesis; 

ui  "  t  “  ~ct  |  %  (  5  ) 

where,  following  Johnson  [111,  the  constant  c  was  taken  to 
be  0.35.  This  relatively  simple  formulation  is  at  least  a 
more  generally  valid  statement  than  eq(3)  since  it  gives  rise 
to  non-isotropic  diffusivities  and,  indeed,  can  produce  heat 
fluxes  in  directions  where  there  is  no  temperature  variation, 
j A  recent  experimental  paper  by  Tavoularis  and  Corrsin  [12] 
has  helped  underline  the  importance  of  this  feature. 

! 

The  turbulent  scalars  k  and  e  are  obtained  from  the 
I  usual  transpoi^  equations  which  in  toroidal  coordinates 
become;  pC(k)  .  D(k)  +  sD(k)  +  P  _  e  (6) 

.  pC(e)  -  D(e)  +  SD(e)  +  cel  ^  “  ce2  7"  ( 7  > 


where  the  C  operator  is  the  same  as  in  eq(l)  while,  in  the 
case  of  the  ASM  treatment,  the  D  and  S_  operators  take  a 
cannon  form  for  both  k  and  e:  u 
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The  diffusion  coefficients  for  k  and  e  (i.e.  c.  and  cE)  are 
assigned  the  values  0.22  and  0.15,  again  in  line  with 
standard  practice. 


I  In  the  case  of  the  EVM  the  diffusion  'sources’  of  k  and 
;  e*  SD(v),  are  2ero  and  the  diffusion  operator  D(*)  takes  the 
! same  form  as  for  the  mean  velocity  component  with  u  replaced 
by  (  v  +  pvfc/ ).  The  turbulent  Prandtl  numbers  for 
,  turbulent  energy  ( a  )  and  energy  dissipation  (®E)  take  the 
j values  1.0  and  1.22. 

|  The  turbulent  energy  generation  rate  P  is  given  by: 

I 
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When  the  EVM  is  adopted,  eq(2)  is  substituted  into  (8)  so 
,  that  only  mean  velocity  gradients  and  the  turbulent  viscosity 
arise  in  evaluating  P. 

;  (b)  The  low-Reynolds-number  region  :  In  the  immediate 
vicinity  of  the  wall  viscous  effects  on  the  turbulent 
transport  processes  become  substantial.  Although  low- 
:  Reynolds- nisnber  forms  of  the  k  ~  c  EVM  were  available  (e.g. 
i  [13] } ,  the  handling  of  the  k  and  e  equations  in  the  viscosity 
j  affected  layers  is  fairly  tricky.  Thus,  in  view  of  the  heavy 
burden  that  our  calculations  would  inevitably  place  on  our 
!  computing  resources,  a  simpler  scheme  was  adopted.  In  axi- 
symmetric  boundary  layer  flows  near  spinning  surfaces  it  is 
known  [14]  that  Prandtl' s  mixing  length  hypothesis  with  Van 
Driest's  [15]  damping  function  leads  to  a  broadly 
satisfactory  description  of  the  momentum  and  heat  transport 
processes.  Accordingly,  in  the  annular  ring  0.96  <  r /a  4  1.0 
it  is  assumed  that 
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where  the  mixing  length  f  varies  with  distance  y  from  the 
pipe  wall  as 


im  *  0.419  y  (1  -  exp  -  (p  yUT/  26  lO) 


The  turbulent  Prandtl  number  for  heat  transfer  is  again  taken 
as  0.9. 


3.  DISCRETIZATION  AND  SOLUTION 
3.1  General  Strategy 

A  complete  presentation  of  the  procedure  developed  for 
solving  the  system  of  equations  introduced  in  §2  is  given  by 
Iacovides  [16] .  In  the  next  sub-section  we  focus  on  the 
methods  for  handling  the  Reynolds  stresses  when  the  ASM  is 
adopted;  these  represent  the  principal  innovations  of  this 
paper  so  far  as  the  numerics  are  concerned.  Although  the 
introduction  of  the  ASM  has  necessitated  a  major  re¬ 
structuring  of  the  code,  the  same  broad  solution  strategy 
evolved  for  use  with  EVMs  (and  described  in  more  detail  in 
refs[4],  [5]  and  [6])  has  been  retained.  Namely; 

•  A  finite-volume,  primitive-variable  decomposition 
using  the  now  customary  staggered  mesh  in  which  the 
nodes  for  velocity  are  located  on  the  boundaries  of 
the  scalar  nodes  (p,  T,  k  and  e ) . 

•  The  use  of  QUICK  [17]  for  approximating  convective 
transport  of  all  variables  in  the  cross-stream 
(r~$)  plane. 

•  The  adoption  of  a  'semi-elliptic'  treatment  in 
which  all  dependent  variables  except  p  are  stored 
on  just  two  adjacent  r—  4>  planes,  with  values  being 
successively  transferred  from  the  downstream  to  the 
upstream  register  as  each  streamwise  iteration  is 
made. 

•  The  use  of  Patankar's  SIMPLER  algorithm  [18]  for 
obtaining  the  pressure  field  in  compliance  with  the 
continuity  equation.  Some  adaptation  of  the 
original  scheme  was  needed  to  incorporate  it  in  a 
semi-elliptic  framework. 


The  employment  of  PSL  [19]  over  a  thin  band  near 
the  wall  (which  for  convenience  we  take  to  coincide 
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with  the  region  where  the  mixing-length  scheme  is 
used  for  computing  turbulent  viscosities)  in  which 
the  pressure  is  obtained  from  radial  equilibrium 
(rather  than  continuity)  and  where,  in  consequence, 
three-dimensional  storage  of  p  is  not  required. 

The  180°  U-bend  computations  (both  ASM  and  EVM)  have 
evaluated  the  dependent  variables  on  a  24  (azimuthal)  x  28 
(radial)  x  160  (streamwise)  mesh.  Of  the  radial  nodes,  8 
were  located  in  the  PSL  region.  Convergence  was  flagged  when 
the  sum  of  the  absolute  mass  residuals  for  each  control 
volume  was  reduced  below  1%.  This  degree  of  convergence  was 
secured  after  approximately  100  streamwise  iterations  over 
the  domain. 

3.2  Treatment  of  the  Reynolds  Stresses 

(a)  Location  of  the  stresses  :  The  Reynolds  stresses  lie  on 
a  mesh  staggered  with  respect  to  the  velocity  nodes  as 
indicated  in  Figure  2.  This  choice  promotes  stability  and 
keeps  to  a  minimum  the  amount  of  interpolation.  Notice,  for 
example,  that,  located  where  it  is,  the  shear  stress  vw  lies 
on  the  boundary  of  both  a  W  control  volume  and  a  V  control 
volume  where  it  is  needed  in  applying  the  force  momentum 
principle  to  these  velocity  components.  Likewise  vw  is 
generated  (in  part)  by  the  differences  in  both  the  V  and  W 
velocities  at  the  nodes  straddling  the  shear  stress. 

(b)  Pseudo-viscosity  tensor  and  residual  stress  :  Over  the 
whole  region  In  which  the  ASM  is  applied  the  viscous 
diffusion  terms  play  an  insignificant  roles  the  equations  are 
source-term  dominated,  the  sources  (and  sinks)  being  the 
Reynolds  stresses.  In  an  EVM  treatment,  of  course,  the 
stresses  are  re-expressed  in  terms  of  a  viscosity  times  a 
mean  strain,  introducing  into  the  difference  equations  a 
strong  coupling  between  the  velocities  at  one  point  and  those 
at  surrounding  nodes.  7his  linkage  decisively  improves  the 
stability  of  the  equation  set. 

With  an  ASM  the  intercoupling  among  the  sources  is 
vastly  more  complex  and,  based  on  experience  in  much  simpler 
two-dimensional  flows,  attempts  to  solve  the  equations  by 
treating  the  Reynolds  stresses  entirely  as  explicit  sources 
and  sinks  would  not  have  led  to  convergence.  Accordingly,  a 
(fairly)  strong  implicit  coupling  among  neighbouring 
velocities  is  created  by  introducing  a  pseudo-viscosity 
tensor  In  the  U,  momentum  equation,  the  value  of  viy  is 

simply  *Jie  magnitude  of  the  coefficient  multiplying  MJj/axj 
in  the  ASM  expression  for  UjU..  For  the  toroidal  coordinate 
system  adopted,  the  components  of  the  pseudo-viscosity  are 
given  in  Table  2.  The  part  of  the  Reynolds  stress  not 
accounted  for  through  the  pseudo-viscosity  is  termed  the 
"residual  stress",  uju^,  the  components-  of  which,  for 


Table  3  Residual  Stress  Tensor  u^u 
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completeness,  are  given  in  Table  3.  The  residual  stresses 
can  be  accommodated  as  sources  and  sinks  in  the  equations 
without  provoking  divergence. 


(c)  Generation  term  treatment  :  Although  it  is  physically 
possible  tor  interactions  among  the  turbulent  stresses  and 
mean  strains  to  lead  to  negative  "generation"  rates  of 
turbulence  energy,  in  practical  shear  flows  this  occurs  only, 
if  at  all,  in  limited  regions  and  the  magnitude  of  this  extra 
sink  of  turbulence  energy  is  so  small  as  to  be  an  unimportant 
contributor  to  the  energy  balance  in  the  region  in  question. 
In  an  iterative  numerical  solution  there  is  no  guarantee  that 
during  the  progress  of  iteratidn  substantial  negative 
generation  rates  will  not  occur,  however.  When  they  do,  a 
very  rapid  divergence  can  result.  This  behaviour  has  been 
reported  by  others  who  have  usually  restored  convergence  by 
setting  P  to  zero  whenever  negative  values  arise.  This 
practice  proved'  ineffective  in  the  present  study,  however . 
The  problem  arises  not  just  from  the  kinetic  energy  and 
dissipation  rate  equations  where  P  appears  as  a  source.  The 
anisotropy  of  the  Reynolds  stresses  is  directly  sensitive  to 
the  value  of  P/e  and  it  is  this  sensitivity  that  is  the  most 
damaging  to  securing  convergence .  The  problem  has  been 
removed  by  the  following  treatment.  We  define  a  residual 
turbulence  energy  generation  P  as 


P  =  P  _  PEVM 

where  P__.  is  the  value  of  the  turbulence  energy  generation 
given  cyan  eddy-viscosity  approximation  of  the  Reynolds 
stresses,  viz 
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While  Pp^  is  always  positive,  "p  may  be  positive  or  negative. 
In  evaluating  the  Reynolds  stresses  from  the  ASM  equations 
the  dimensionless  parameter  P/  e  is  replaced  by  P’/e'  where 

p'  -  W  Hax  (5'0) 

e  *  -  e  +  Max  (-5,0) 

Thus,  when  P  is  positive,  P'/C  is  identical  to  P/e  while 
negative  residual  generations  augment  the  dissipation  term. 
In  the  kinetic  energy  equation  P'  and  C  likewise  replace  P 
and  e  while,  in  the  dissipation  rate  equation,  if  P  is 
negative  the  source/sink  terms  become 
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Although  algebraically  the  k  and  e  equations  are  unchanged  by 
this  rearrangement,  there  are  implied  differences  in  the 
numerical  processing  since,  following  long  established 
practice,  source  terms  are  handled  explicitly  while  sink 
terms  are  treated  implicitly. 

(d)  Solution  of  the  algebraic  stress  equations  :  Perhaps  the 
obvious  way  to  handle  die  six  coupled  algebraic  equations  for 
the  Reynolds  stresses  would  be  to  solve  the  6x6  matrix  at 
each  point  in  the  flaw.  One  may  recognize,  however,  that 
there  are  distinct  differences  between  the  role  of  the 
diagonal  and  off-diagonal  elements.  For  a  given  node 
location  the  three  diagonal  elements  (the  normal  stresses) 
are  located  at  a  coincident  point  while  the  off-diagonal 
terms  (the  shear  stresses)  are  scattered.  Moreover,  the 
normal  stresses  can  physically  never  be  negative  -  and  must 
not  be  allowed  to  become  so  during  progress  towards  a 
solution  for  that  would  produce  negative  components  of 
pseudo-viscosity  (Table  2) . 

For  these  and  other  reasons  we  have  followed  the 
practice  of  Huang  and  Leschziner  (20]  in  handling  the  normal 
and  shear  stresses  separately < T\  However,  because  of  the 
strong  coupling  between  the  shear  and  normal  stresses,  direct 
solution  of  the  two  3x3  matrices  is  not  viable.  Instead,  a 
successive-substitution  approach  has  been  adopted.  lacovides 
[16]  provides  a  detailed  account  of  the  steps;  here  we  simply 
sirrmarize  the  main  points.  At  every  point  on  an  c~ 4>  plane 
the  three  normal  stresses  ace  obtained,  explicitly  using 
values  of  (the  other)  Reynolds  stresses  from  the 
corresponding  node  on  the  upstream  plane  and  current-plane 
values  of  other  dependent  variables.  Initially  just  that 
single  iteration  is  made  at  any  plane;  this  limitation  is 
applied  not  just  because  there  is  little  merit  in  obtaining 
very  precise  values  for  the  turbulent  stresses  while  the 
momentum  and  continuity  equations  are  far  from  converged,  but 
also  because  such  attempts  to  obtain  more  accurate  stresses 
actually  slow  convergence  -  and  on  occasions  have  caused 
divergence.  'Ihe  three  shear  stresses  are  determined  next  in 
much  the  same  way.  As,  with  successive  streamwise  passes 
over  the  solution  domain,  the  residuals  in  the  mean  field 
variables  become  small,  first  one  and  then  two  additional 
iterations  are  made  on  the  Reynolds  stresses  at  each  point. 
Because  of  the  large  number  of  streamwise  planes  (160)  the 
change  in  Reynolds  stress  over  a  single  forward  step  is 
usually  only  a  few  percent.  Thus,  the  use  of  the  upstream 
value  gives  an  initial  estimate  of  each  stress  component  that 
is  quite  close  to  the  actual  solution;  only  a  few  iterations 


I.  Huang  and  Leschziner  have  been  concerned  with  two- 
dimensional  flows  in  wiiich  only  one  shear  stress 
component  is  non-zero. 
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therefore  required.  "  '  ~  | 
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3.3  Boundary  Conditions 

I 

For  the  case  of  the  U-bend,  the  velocity  and  turbulence 
profiles  at  inlet  appropriate  to  fully-developed  flow  were 
prescribed.  These  were  obtained  by  running  the  code  for  tne 
case  of  a  straight  pipe.  At  exit  the  only  constraint  applied 
was  on  pressure  where  the  streamwise  gradient  was  taken 
uniform  over  the  pipe  cross  section. 

At  the  pipe  wall  all  velocity  components  were  set  to 
zero.  The  turbulence  energy  and  its  dissipation  rate  were 
not  solved  in  the  low-Reynolds-number  region.  However,  at 
the  interface  between  this  sublayer  and  the  fully  turbulent 
core,  effective  values  of  k  and  e  were  obtained  by  setting 
equal  to  the  turbulence  energy  generation  rate  (given  by 
the  mixing  length  hypothesis)  and  requiring  that  the 
viscosity  given  by  the  mixing  length  and  formulae  should 
be  the  same.  -  i 
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The  value  of  the  circunferential  velocity  U  was  set  to 
zero  along  the  syninetry  plane  while  the  gradient  of  all  other 
dependent  variables  was  set  to  zero.  References  [4]  and  [5] 
describe  the  procedure  for  handling  the  symmetry  at  the  pipe 
centre. 

4.  PRESENTATION  AND  DISCUSSION  OF  RESULTS 

Comparisons  are  here  limited  to  the  180°-bend 
experiments  undertaken  by  Azzola  and  Humphrey  [21]  and  the 
authors  [22]  and  to  the  heat  transfer  data  of  flow  in  a  coil 
by  Seban  and  McLaughlin  [23].  A  fuller  set  of  comparisons 
is  given  by  Iacovides  [16].  The  reader  is  referred  to  ref 
[4]  and  [5]  which  show  excellent  agreement  with  data  of 
various  laminar  flows,  thus  providing  an  indication  that 
numerical  errors  are  small  with  the  discretization  practices 
and  mesh  densities  here  adopted.  Although  the  turbulent  flow 
transport  equations  using  an  ASM  are  far  more  complicated 
than  for  laminar  flow,  substantial  grid  refinement  for  the 
Seban-McLaughlin  test  case  (see  below)  produced  insignificant 
changes  in  heat  transfer  coefficient. 

The  Reynolds  numbers  of  the  U-bend  computations  do  not 
correspond  exactly  with  those  of  the  experiments  principally 
because  computational  work  has  proceeded  in  parallel  with 
experiments  anc.  'nominal'  flow  rates  were  adopted.  At  a 
Reynolds  number  of  6  x  10 4  a  difference  of  10%  in  Re  is  known 
not  to  have  a  perceptible  effect  on  the  flow  structure  (ID, 


II.  Azzola  and  Humphrey  [6]  report  measurements  for  two 
values  of  Re;  57,400  and  110,000.  Differences  in 
velocity  profiles  at  these  two  Reynolds  numbers  are  in¬ 
distinguishable  from  experimental  uncertainty. 


Fig.  4  ASM  and  EVM  streamwise  velocity  profiles  along 
surfaces  parallel  to  plane  of  symmetry 

a)  8  =  90°  b)  8  =  135° 


Fig.  6  ASM  vector  plots  of  secondary  velocity.  Arrow 

length  such  that  10%  of  pipe  diameter^  corresponds 
to  16%  of  bulk  streamwise  velocity,  W, 

a)  0  =  45°  b)  135° 
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Predicted  levels  of  Nusselt  number  have  been  converted  to  the 
'true*  Reynolds  number  by  multiplying  by  the  ratio  of  the 
true  to  the  computed  Reynolds  number  to  the  power  O.o. 

Azzola  and  Humphrey  [21]  report  velocity  profiles  along 
the  radial  line  for  $  =  90°,i.e.  orthogonal  to  the  plane  of 
symmetry.  Their  measurements  of  the  streamwise  velocity  to 
are  compared  in  Figure  3  with  predictions  of  both  the  EVM  and 
ASM  schemes.  Both  models  capture  the  overshoot  of  the  mean 
velocity  by  90°  into  the  bend  and  the  subsequent  flattening 
of  the  velocity  profile  over  the  central  core.  As  it  turns 
out,  the  velocity  distribution  along  this  line  is  not  a 
sensitive  discriminator  between  the  two  models.  Figure  4 
therefore  shows  computed  mean  velocity  distributions  along 
lines  parallel  to  the  plane  of  symmetry  of  the  flow  at  90° 
and  135°  into  the  U-bend.  The  ASM  distributions  show  a  more 
distinct  secondary  peak  near  the  inside  of  the  bend  titan 
those  generated  by  the  EVM.  Based  on  our  experience  of  flow 
in  90°  bends  [6J  we  have  no  doubt  that  the  sharper  peaking 
exhibited  by  the  ASM  computations  better  accords  with  the 
experimental  behaviour. 

The  secondary  ^circumferential)  velocity  profile  along 
the  90°  line  is  compared  with  experiment  at  several  stations 
in  Figure  5.  Both  predicted  curves  indicate,  in  accord  witn 
experiment,  £hat  the  peak  secondary  flow  occurs  at  45°  and 
that  by  135°  the  flow  near  the  pipe  centre  is  directed 
towards  the  inside  of  the  bend.  The  overall  implications  of 
these  changes  in  the  profile  are  most  clearly  seen  by  the 
secondary  flow  vector  plots  of  Figure  6:  the  conventional 
picture  of  £he  single-cell  secondary  circulation  at  45°  gives 
way  at  135°  to  a  kidney-shaped  primary  eddy  and  a  weax 
counter-rotating  secondary  vortex  near  the  tube  centre.  The 
ASM  predictions  in  Figure  5  achieve  somewhat  closer  agreement 
with  the  measured  profiles  than  do  the  EVM  calculations.  The 
largest  relative  differences  between  the  two  numerical 
computations  actually  occur  downstream  of  the  bend,  the  ASM 
results  showing  a  stronger  movement  towards  (what  had  been) 
the  inside  of  the  bend  than  near  the  wall,  the  reverse  of 
what  the  EVM  predicts. 

Turning  now  to  the  convective  heat  transfer  results,  the 
Nusselt  number  profiles  for  the  Seban-McLaughlin  coil  shown 
in  Figure  7  display  a  2:1  variation  in  heat  transfer 
coefficient  between  the  inside  and  outside  of  the  coil 
arising  partly  from  the  damping/augmentation  of  turbulence 
near  the  convex/concave  surfaces  and  partly  from  the 
secondary  flow.  Since,  for  this  case,  the  flow  is  fully 
developed  in  the  streamwise  direction,  storage  of  the 
pressure  field  imposes  no  large  core  requirements. 
Accordingly,  a  second  ASM  computation  was  made  with  40  radial 
and  30  circumferential  nodes.  The  calculated  levels  of  Nu 
differed  by  at  most  1*>%  from  those  obtained  with  the  standard 


28  x  2 4  mesh.  The  result  helps  confirm  our  view  that 
numerical  errors  have  been  kept  to  unimportant  levels.  The 
circumferential  variation  of  Nusselt  number  computed  with  tne 
ASM  is  somewhat  less  than  the  measured  though  agreement  is 
better  than  for  the  EVM  results. 

In  fact,  this  set  of  measurements  do  not  provide  an 
ideal  test  case  because  experimental  difficulties  give  rise 
to  uncertainties  in  the  data  comparable  with  the  differences 
shown  between  experiment  and  computation  in  Figure  7.  The 
principal  convective  heat  transfer  test  is  tnus  the  data  of 
developing  'flow  in  a  180°  bend  obtained  recently  at  UMIST 
[22].  The  experiment  is  one  at  uniform  wall  temperature,  the 
local  heat  flux  being  obtained  by  measuring  the  electrical 
power  supply  to  a  small  circular  disc  required  to  keep  the 
disc  at  the  same  temperature  as  the  rest  of  the  thick 
aluminium  U-bend.  Heating  starts  30  diameters  upstream  of 
the  bend  so  the  thermal  field  is  nearly  fully  developed  at 
the  start  of  the  bend.  Hie  experimental  uncertainty  of  the 
resultant  Nusselt  numbers  is  believed  to  be  within  ±7%.  The 
geometric  arrangement  of  the  apparatus  is  a  (roughly  twice) 
scale  replica  of  the  U-bend  of  Azzola  and  Humphrey  and  the 
test  flow  Reynolds  numbers  in  the  two  cases  differ  by  no  more 
than  5%.  Moreover,  the  difference  in  temperature  between  the 
wall  and  the  entering  fluid  was  only  20°C  so  the  thermal 
experiment  was  essentially  one  with  uniform  fluid  properties. 
We  may  thus  assume  that  the  velocity  field  measured  by  Azzola 
and  Humphrey  [21]  describes  the  hydrodynamic  development  in 
the  thermal  experiment  also. 

Normalized  temperature  profiles  measured  with  a  thermo¬ 
couple  rake  are  shown  in  Figure  8  for  three  positions  in  the 
bend  and  at  6  diameters  downstream.  The  secondary  circulation 
carries  hot  fluid  from  near  the  inside  of  the  bend  along  the 
symmetry  plane  (viz.  Figure  6)  and  this  leads  to  an  ‘over¬ 
shoot’  in  the  temperature  profile  The  pattern  is  well 
reproduced  by  the  two  sets  of  predictions  though,  in  so  far 
as  differences  are  discernible,  the  ASM  computations  give  the 
better  agreement  with  the  measurements. 

The  development  of  the  local  Nusselt  number  around  the 
bend  is  compared  in  Figure  9  with  predictions  generated  by 
the  two  turbulence  models;  here,  because  the  experiment  was 
one  of  constant  wall  temperature,  the  variation  of  the  bulk 
mean  temperature  around  the  bend  is  unknown.  Tne  Nusselt 
number  presented  in  Figure  9  is  thus  based  or.  the  difference 
between  the  wall  temperature  and  that  at  the  centre  of  the 
pipe.  Hie  fully-developed  flow  Nusselt  number  for  a  straight 


III. Hie  normalization  makes  the  coolest  fluid  correspond  with 
the  maximum  in  the  temperature  profile. 


Fig.  9  Development  of  local  Nusselt  number  around  180  bend 
Key  as  Fig.  8 


tube  at  this  Reynolds  number  is  about  110.  As  the  flow 
develops  around  the  bend  the  Nusselt  number  is  depressed  on 
the  inside  of  the  bend  and  augmented  on  the  outside.  'Me 
maximum  level  occurs  at  the  90°  station  with  a  ratio  of 
maximum sminimum  Nu  of  about  3.5:1.  Over  the  second  half  of 
the  bend  there  is  a  gradual  reduction  in  this  ratio  and  in 
the  average  level.  At  55°  the  two  predicted  curves  show 
nearly  the  same  distribution  reproducing  fairly  well  the 
measurements  through  with  levels  about  10%  lower.  By  90 
significant  differences  have  developed  both  between  the 
predictions  themselves  and  the  measurements  and  predictions. 
Me  Nu  levels  given  by  the  EVM  are  substantially  too  low  over 
the  whole  circumference  -  nearly  30%  low  near  the  outside  of 
the  bend.  Me  ASM  results  accord  more  closely  with  the  data, 
particularly  over  the  inside  of  the  bend,  but  they  still  give 
values  of  Nu  nearly  25%  below  the  measurements  at  t£e 
outside.  A  similar  pattern  is  exhibited  at  125  and  165  . 
At  the  former  position  we  note  that  the  EVM  predicts  a 
minimum  Nu  off  the  symmetry  plane  but  the  ASM  curve  is  in 
markedly  better  agreement  with  experiment,  the  maximum  dis¬ 
agreement  with  experiment  now  being  only  15%.  Downstream  of 
the  bend  the  ASM  predictions  display  a  faster  decay  towards 
fully  developed  flow  than  the  EVM  and  in  Figure  9(c)  this 
model  predicts  levels  of  Nu  slightly  below  the  EVM  values. 
We  have  seen  earlier  that  the  secondary  flow  behaviour 
(Figure  5)  indicated  that  the  ASM  computations  predict  more 
correctly  the  downstream  recovery  so  the  'good'  performance 
of  the  EVM  at  this  station  is  really  a  fortuitous  outcome  of 
a  weakness  in  the  model. 

Finally,  we  must  draw  attention  to  an  interesting 
anomaly  in  the  predictedQheat  transfer  behaviour  at  the  entry 
to  the  bend.  At  15°  into  the  bend.  Figure  10,  the 
experimental  data  indicate  that  the  heat  transfer  coefficient 
on  the  inner  line  of  symmetry  is  more  than  40%  below  that  on 
the  outside  whereas  the  predictions  show  an  essentially 
uniform  level.  Mis  might  seem  to  suggest  that  the 
computations  were  not  responding  rapidly  enough  to  the 
effects  of  applying  curvature  of  the  streamlines.  However, 
the  streamwise  acceleration  on  the  inside  of  the  bend  is 
rather  severe  at  the  bend  entry:  tha-parameter  -  (v/u^)3P/r3e 
reaches  a  value  in  excess  of  2  x  10~  .  If  this  strength  of 
acceleration  were  applied  to  a  two-dimensional  boundary  layer 
laminarization  would  occur.  When  applied  for  only  a  short 
distance  complete  reversion  to  laminar  flow  does  not  occur 
but  the  viscous  sublayer  becomes  substantially  thicker  than 
normal  substantially  increasing  thermal  resistance  and  thus 
reducing  the  heat  transfer  coefficient.  On  the  outside  of 
the  bend  there  is  a  substantial  adverse  pressure  gradient 
that  will  tend  to  thin  the  sublayer  producing  the  reverse 
effect  to  that  found  on  the  inside  of  the  bend.  Not 
unexpectedly  our  computations  based  on  the  standard  Van 
Driest  mixing  length  treatment  for  the  near-wall  region 


entirely  miss  this  unforeseen  (end  hitherto  unreported) 
phenomenon  at  entry  to  the  bend.  It  is,  of  course,  possible 
to  tune  the  mixing  length  hypothesis  to  Reproduce  at  least 
partially  this  laminar ization  by  making  A  a  function  of  the 
local  acceleration  (this  has  been  done  in  many  nuxing-lengch 
variants  of  the  late  '60s).  A  better  approach,  nowever, 
would  probably  be  to  employ  a  low-Reynolds-number  version  of 
the  2-equation  EVM,  for  example  [13],  which  is  known  to  do 
fairly  well  at  predicting  laminarization  without  introducing 
flow-dependent  parameters  into  the  model  coefficients. 

5.  CONCLUSIONS 

The  paper  has  reported  the  various  tactics  adopted  to 
achieve  convergent  numerical  predictions  of  strongly  three- 
dimensional  flow  incorporating  an  algebraic  stress  model  of 
turbulence  and  has  presented  predictions  for  flow  in  a  180° 
U-bend.  Apart  from  the  partial  laminarization  that  appears 
to  occur  at  entry  to  the  bend,  overall,  the  ASM-based 
computations  have  achieved  a  satisfactory  simulation  of  the 
resultant  flow  and  thermal  fields.  Indeed,  this  model 
achieves  consistently  better  results  than  the  EVM,  albeit  at 
the  expense  of  roughly  double  the  computational  effort.  A 
remaining  weakness  of  the  model  is  that  heat  transfer 
coefficients  are  predicted  significantly  too  low  on  the 
outside  (concave  surface)  of  the  bend.  A  more  detailed 
experimental  mapping  of  the  mean  and  turbulent  velocity  field 
in  this  flow  would  be  very  helpful  in  tracing  the  source  of 
this  disagreement  in  the  model. 
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